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Abstract We consider time-periodically perturbed ID Hamiltonian systems pos- 
sessing one or more separatrices. If the perturbation is weak, then the separatrix 
chaos is most developed when the perturbation frequency lies in the logarithmically 
small or moderate ranges: this corresponds to the involvement of resonance dynam- 
ics into the separatrix chaos. We develop a method matching the discrete chaotic 
dynamics of the separatrix map and the continuous regular dynamics of the reso- 
nance Hamiltonian. The method has allowed us to solve the long-standing problem 
of an accurate description of the maximum of the separatrix chaotic layer width as 
a function of the perturbation frequency. It has also allowed us to predict and de- 
scribe new phenomena including, in particular: (i) a drastic facilitation of the onset 
of global chaos between neighbouring separatrices, and (ii) a huge increase in the 
size of the low-dimensional stochastic web. 
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1 Introduction 



Separatrix chaos is the germ of Hamiltonian chaos f5T\. Consider an integrable 
Hamihonian system possessing a saddle, i.e. a hyperboHc point in the one-dimensional 
case, or a hyperbolic invariant torus, in higher-dimensional cases. The stable (in- 
coming) and unstable (outgoing) manifolds of the saddle are called separatrices 
ifTSl : they separate trajectories that have different phase space topologies. If a weak 
time-periodic perturbation is added, then the separatrix is destroyed; it is replaced 
by a separatrix chaotic layer (SCL) ||5Tl[T8ll23ll29l . Even if the unperturbed system 
does not possess a separatrix, the resonant part of the perturbation generates a sepa- 
ratrix in the auxiliary resonance phase space while the non-resonant part of the per- 
turbation destroys this separatrix, replacing it with a chaotic layer ISTl [TSl l23l [TOl . 
Thus separatrix chaos is of a fundamental importance for Hamiltonian chaos. 

One of the most important characteristics of SCL is its width in energy (or ex- 
pressed in related quantities). It can be easily found numerically by integration of 
the Hamiltonian equations with a set of initial conditions in the vicinity of the sep- 
aratrix: the space occupied by the chaotic trajectory in the Poincare section has a 
higher dimension than that for a regular trajectory, e.g. in the 3/2D case the regular 
trajectories lie on lines i.e. ID objects while the chaotic trajectory lies within the 
SCL i.e. the object outer boundaries of which limit a 2D area. 

On the other hand, it is important to be able to describe theoretically both the 
outer boundaries of the SCL and its width. There is a long and rich history of the 
such studies. The results may be classified as follows. 



1.1 Heuristic results 

Consider a ID Hamiltonian system perturbed by a weak time-periodic perturbation: 



where Ho{p,q) possesses a separatrix and, for the sake of notational compactness, 
all relevant parameters of Hq and V, except possibly for cOf, are assumed to be ~ 1 . 

Physicists proposed a number of different heuristic criteria lf53llTb. 23 , 55 . 5111521 
for the SCL width AE in terms of energy E = HQ{p,q) which gave qualitatively 
similar results: 



H = Ho{p,q)+hV{p,q,t), 
V{p,q,t + 27z/(Of)^V{p,q,t), 



/l< 1 



(1) 



AE = AE{cOf) - (0/5, 
d = h\e\, 



(2) 




|£| oc exp{-a(Of) < 1 



(fl 



1) 



for CO/ ~ 1 , 
for 0)/ > L 
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The quantity 5 = /2|e| is called the separatrix split 1511 (see also Eq. (4) below): 
it determines the maximum distance between the perturbed incoming and outgoing 
separatrices I53[l0l|23]|55l|5ll|52l[ll[l8l|29l. 

It follows from (2) that the maximum of AE should he in the frequency range 
(Of ^ I while the maximum itself should he ^ h: 

/i£max = max{AE{(Of)} ~ h, cof""""^ ~ 1. (3) 



1.2 Mathematical and accurate physical results 

Many papers studied the SCL by mathematical or accurate physical methods. 

For the range cOf^l, many works studied the separatrix splitting (see the review 
IfTSlI and references therein) and the SCL width in terms of normal coordinates (see 
the review 1291 and references therein). Though quantities studied in these works 
differ from those typically studied by physicists l53l[T0ll23ll55ll5Tll52l . they im- 
plicitly confirm the main qualitative conclusion from the heuristic formula (2) in the 
high frequency range: provided that O/ S> 1 the SCL width is exponentially small. 

There were also several works studying the SCL in the opposite (i.e. adiabatic) 
limit (Of 0: see e.g. l27l fT4l l28l l42l l45l and references therein. In the context 
of the SCL width, it is most important that AE{cOf Q) ^ h for most of the sys- 
tems l27l [T4l l28l . For a particular class of systems, namely for ac-driven spatially 
periodic systems (e.g. the ac-driven pendulum), the width of the SCL part above 
the separatrix diverges in the adiabatic limit l42l |45| : the divergence develops for 
co/<l/ln(l//i). 

Finally, there is a qualitative estimation of the SCL width for the range (Of ^ I 
within the Kolmogorov-Arnold-Moser (KAM) theory l29l . The quantitative esti- 
mate within the KAM theory is lacking, apparently being very difficult for this fre- 
quency range |fT7|. It follows from the results in l29\ that the width in this range is 
of the order of the separatrix split, which itself is of the order of h. 

Thus it could seem to follow that, for all systems except ac-driven spatially peri- 
odic systems, the maximum in the SCL width is ~ h and occurs in the range cOf ^ I, 
very much in agreement with the heuristic result (3). Even for ac-driven spatially 
periodic systems, this conclusion could seem to apply to the width of the SCL part 
below the separatrix over the whole frequency range, and to the width of the SCL 
part above the separatrix for (Of ^ l/ln(l//i). 
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1.3 Numerical evidence for high peaks in AE{(Of) and their rough 
estimation 

The above conclusion disagrees with several numerical studies carried out during 
the last decade (see e.g. EH |45l [3l|25l |40l |24l gTl ED) which have revealed the 
existence of sharp peaks in AE((Of) in the frequency range l/ln(l /h) ^ cOf ^ I the 
heights of which greatly exceed li (see also Figs. 2, 3, 5, 6 below). Thus, the peaks 
represent the general dominant feature of the function AE{(Of). They were related 
by the authors of 1341 125] |40] |24l |47] |35l to the absorption of nonlinear resonances 
by the SCL. For some partial case, rough heuristic estimates for the position and 
magnitude of the peaks were made in ||34ll35l . 



1.4 Accurate description of the peaks and of the related 
phenomena 

Until recently, accurate analytic estimates for the peaks were lacking. It is explic- 
itly stated in 1241 that the search for the mechanism through which resonances are 
involved in separatrix chaos, and for an accurate analytic description of the peaks 
in the SCL width as function of the perturbation frequency, are being among the 
most important and challenging tasks in the study of separatrix chaos. The first 
step towards accomplishing them was taken through the proposal B3l 1441 of a new 
approach to the theoretical treatment of the separatrix chaos in the relevant fre- 
quency range. It was developed and applied to the onset of global chaos between 
two close separatrices. The application of the approach l43l l44ll to the commoner 
single-separatrix case was also discussed. The approach has been further developed 
l38l[39l . including an explicit theory for the single-separatrix case ll39l . 

The present paper reviews the new approach l43l l44l [38l [39l and its applications 
to the single-separatrix 1391 and double-separatrix B3l l44l cases. We also briefly 
review application to the enlargement of the low-dimensional stochastic web l46l 
and discuss other promising applications. 

Though the form of our treatment differs from typical forms of mathematical the- 
orems in this subject (cf. |fT8ll29l ). it yields the exact expressions for the leading term 
in the relevant asymptotic expansions (the parameter of smallness is a = 1 / ln( 1 //j)) 
and, in some case, even for the next-order term. Our theory is in excellent agreement 
with results obtained by numerical integration of the equations of motion. 

Sec. 2 describes the basic ideas underlying the approach. Sec. 3 is devoted to the 
leading-order asymptotic description of the single-separatrix chaotic layers. Sec. 4 
presents an asymptotic description of the onset of global chaos in between two close 
separatrices. Sec. 5 describes the increase in sizes of a stochastic web. Conclusions 
are drawn in Sec. 6. Sec. 7 presents the Appendix, which explicitly matches the 
separatrix map and the resonance Hamiltonian descriptions for the double-separatrix 
case. 
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2 Basic ideas of the approach 

The new approach Il43ll44l[38l[39l may be formulated briefly as a matching between 
the discrete chaotic dynamics of the separatrix map in the immediate vicinity of the 
separatrix and the continuous regular dynamics of the resonance Hamiltonian be- 
yond that region. The present section describes the general features of the approach 
in more detail. 

Motion near the separatrix may be approximated by the separatrix map (SM) 
ll53][10l|23l|55]|ni|52E|29l[33[35]|4^ This was introduced in 121 and 

its various modifications were subsequently used in many studies. It is sometimes 
known as the whisker map. It was re-derived rigorously in ISTI as the leading-order 
approximation of motion near the separatrix in the asymptotic limit h — 0, and 
an estimate of the error was also carried out in 1311 (see also the review 1291 and 
references therein). 

The main ideas that allow one to introduce the SM ll53l [TOl l23l l55l 1511 l52l [Tl l29l 
|43]|44][3T| are as follows. For the sake of simplicity, let us consider a perturbation 
V that does not depend on the momentum: V = V{q,t). A system with energy close 
to the separatrix value spends most of its time in the vicinity of the saddle(s), where 
the velocity q is exponentially small. Differentiating E = Ho{p,q) with respect to 
time and allowing for the equations of motion of the system (1), we can show that 
dE/dt = dV /dqq °^ q. Thus, the perturbation can significantly change the energy 
only when the velocity is not small i.e. during the relatively short intervals while the 
system is away from the saddle(s): these intervals correspond to pulses of velocity 
as a function of time (cf. Fig. 20 in the Appendix below). Consequently, it is pos- 
sible to approximate the continuous Hamiltonian dynamics by a discrete dynamics 
which maps the energy E, the perturbation angle (p = (Oft, and the velocity sign 
(7 = sgn(^), from pulse to pulse. 

The actual form of the SM may vary, depending on the system under study, but 
its features relevant in the present context are similar for all systems. For the sake of 
clarity, consider the explicit case when the separatrix of H(){p,q) possesses a single 
saddle and two symmetric loops while V — qcos{cOft). Then the SM reads l43l (cf. 
Appendix): 

Ei+i 
(Pi+i 

<7i+I 



: Ei + C7,72esin((j!),), 

a)/7r(3-sgn(£',-+i -£,)) 
•^'■^ 2«(£,+,) 
: a,- sgn(£'., -£■,•+!), = 1, 

e = sgn j^ — 



dt 



dHo 
dp 



sm{(Oft), 



dp 



(4) 
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where is the separatrix energy, Co{E) is the frequency of oscillation with energy 
E in the unperturbed case (i.e. for h = 0), f,- is the instant corresponding to the /-th 
turning point in the trajectory q{t) (cf. Fig. 20 in the Appendix below), and A is 
an arbitrary value from the range of time intervals which greatly exceed the char- 
acteristic duration of the velocity pulse while being much smaller than the interval 
between the subsequent pulses ||53[l0l|23l|55l|5l]|52l[I]|29l[3l]. Consider the two 
most general ideas of our approach. 

1. If a trajectory of the SM includes a state with E = Es and an arbitrary (p and d, 
then this trajectory is chaotic. Indeed, the angle (jO of such a state is not correlated 
with the angle of the state at the previous step of the SM, due to the divergence 
of 0)^' (£ ^ iij). Therefore, the angle at the previous step may deviate from a 
multiple of 2n by an arbitrary value. Hence the energy of the state at the previous 
step may deviate from Es by an arbitrary value within the interval [— /ile|,/!|e|]. 
The velocity sign o is not correlated with that at the previous step eitheo Given 
that a regular trajectory of the SM cannot include a step where all three variables 
change random-like, we conclude that such a trajectory must be chaotic. 
Though the above arguments may appear to be obvious, they cannot be consid- 
ered a mathematically rigorous proof, so that the statement about the chaotic 
nature of the SM trajectory which includes any state with E = Es should be con- 
sidered as a conjecture supported by the above arguments and by numerical it- 
eration of the SM. Possibly, a mathematically rigorous proof should involve an 
analysis of the Lyapunov exponents for the SM (cf. fT^) but this appears to be 
a technically difficult problem. We emphasize however that a rigorous proof of 
the conjecture is not crucial for the validity of the main results of the present pa- 
per, namely for the leading terms in the asymptotic expressions describing (i) the 
peaks of the SCL width as a function of the perturbation frequency in the single- 
separatrix case, and (ii) the related quantities for the double-separatrix case. It 
will become obvious from the next item that, to derive the leading term, it is suffi- 
cient to know that the chaotic trajectory does visit areas of the phase space where 
the energy deviates from the separatrix by values of the order of the separatrix 
split 5 = h\e\, which is a widely accepted fact l53[l0l|23l|55l|5ll|52l[I][l3|29l. 

2. It is well known 123 [IHl El EI] Ell [B 111 |29l [51 [35l that, at the 
leading-order approximation, the frequency of eigenoscillation as function of the 
energy near the separatrix is proportional to the reciprocal of the logarithmic 
factor 



bKCOQ 3-sgn(£-£,) 



In 



AH 



\E-Es\ 

\E-Es\<^AH = Es-Es,, 



^ Formally, sgn(£ — is not defined for E = but, if to shift E from for an infinitesemal 
value, sgn(£ — £,) acquires a value equal to either +1 or —1, depending on the sign of the shift. 
Given that Cj+i is proportional to sgn(£',. — while the latter is random-like (as it has been 

shown above), ffj+i is not con'elated with (7, if i = £, ± 0. 
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where E^t is the energy of the stable states. 

Given that the argument of the logarithm is large in the relevant range of E, 
the function (o{E) is nearly constant for a substantial variation of the argu- 
ment. Therefore, as the SM maps the state {Eq = Es,(Po,Gq) onto the state with 
E =E\= Es + OqIie sin(^o), the value of (o{E) for the given sgn((7o£ sin(9o)) is 
nearly the same for most of the angles (except in the vicinity of multiples of 
n). 



Moreover, if the deviation of the SM trajectory from the separatrix increases fur- 
ther, (£){E) remains close to O)/*' provided the deviation is not too large, namely 
if Ind^ - Es\/h) < \n{AH/h). If (Of ~ Gyr^\ then the evolution of the SM (4) 
may be regular-like for a long time until the energy returns to the close vicinity 
of the separatrix, where the trajectory becomes chaotic. Such behavior is espe- 
cially pronounced if the perturbation frequency is close to 0)^^' or (Or ' or to 
one of their multiples of relatively low order: the resonance between the pertur- 
bation and the eigenoscillation gives rise to an accumulation of energy changes 
for many steps of the SM, which results in a deviation of E from Es that greatly 
exceeds the separatrix split h\e\. Consider a state at the boundary of the SCL. 
The deviation of energy of such a state from Es depends on its position at the 
boundary. In turn, the maximum deviation is a function of (Of. The latter func- 
tion possesses the absolute maximum at (Of close to O)/^^ or O),- ' typicall}0, for 
the upper or lower boundary of the SCL respectively. This corresponds to the 
absorption of the, respectively upper and lower, 1 st-order nonlinear resonance by 



The second of these intuitive ideas has been explicitly confirmed ll43l (see Ap- 
pendix): in the relevant range of energies, the separatrix map has been shown to 
reduce to two differential equations which are identical to the equations of motion 
of the auxiliary resonance Hamiltonian describing the resonance dynamics in terms 
of the conventional canonically conjugate slow variables, action / and slow angle 
i/A = n\i/~(Oft where Xj/ is the angle variable IlOl EH |55l ISB |52l [B (see Eq. (16) 
below) and n is the relevant resonance number i.e. the integer closest to the ratio 
/ (±) 

Thus the matching between the discrete chaotic dynamics of the SM and the con- 
tinuous regular-hke dynamics of the resonance Hamiltonian arises in the following 



^ For the SM relating to ac-driven spatially periodic systems, the time during which the SM un- 
dergoes a regular-like evolution above the separatrix diverges in the adiabatic limit (Of — > 1451 , 
and the width of the part of the SM layer above the separatrix diverges too. However, we do not 
consider this case here since it is irrelevant to the main subject of the present paper i.e. to the 
involvement of the resonance dynamics into the separatrix chaotic motion. 



(o[E)^(or\ 
(ol^^ = (0{Es±h), 



sgn((7oesin((j!)o)) = ±1. 



(6) 



the SCL. 
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way ||43l . After the chaotic trajectory of the SM visits any state on the separatrix, 
the system transits in one step of the SM to a given upper or lower curve in the 
I — ij/ plane which has been called Bsll the upper or lower generalized separatrix 
split (GSS) curve0 respectively: 

£=4ss(r)=£.v±5|sin(v/)|, d = h\e\, (7) 

where 5 is the conventional separatrix split ISTl . e is the amplitude of the Melnikov- 
like integral defined in Eq. (4) above (cf. |l53[l0l|23|55]|in|52l[l][ia|29l|3l|43 
[35]|43]|44l), and the angle iff may take any value either from the range [0, n] or from 
the range [;r,2;r0. 

After that, because of the closeness of cOf to the «-th harmonic of Co{E) in the 
relevant range of £0, for a relatively long time the system follows the nonlinear 
resonance (NR) dynamics (see Eq. (16) below), during which the deviation of the 
energy from the separatrix value grows, greatly exceeding 5 for most of the trajec- 
tory. As time passes, iff is moving and, at some point, the growth of the deviation 
changes for the decrease. This decrease lasts until the system hits the GSS curve, 
after which it returns to the separatrix just for one step of the separatrix map. At 
the separatrix, the slow angle iff changes random-like, so that a new stage of evolu- 
tion similar to the one just described occurs, i.e. the nonlinear resonance dynamics 
starting from the GSS curve with a new (random) value of iff. 

Of course, the SM cannot describe the variation of the energy during the velocity 
pulses (i.e. in between instants relevant to the SM): in some cases this variation can 
be comparable to the change within the SM dynamics. This additional variation will 
be taken into account below, where relevant. 

One might argue that, even for the instants relevant to the SM, the SM describes 
the original Hamiltonian dynamics only approximately lISTI and may therefore miss 
some fine details of the motion: for example, the above picture does not include 
small windows of stability on the separatrix itself. However these fine details are 
irrelevant in the present context, in particular the relative portion of the windows of 
stability on the separatrix apparently vanishes in the asymptotic limit /z — > 0. 

The boundary of the SM chaotic layer is formed by those parts of the SM chaotic 
trajectory which deviate from the separatrix more than others. It follows from the 
structure of the chaotic trajectory described above that the upper/lower boundary 
of the SM chaotic layer is formed in one of the two following ways 14311441 : (i) if 
there exists a self-intersecting resonance trajectory (in other words, the resonance 
separatrix) the lower/upper part of which (i.e. the part situated below/above the 



' The GSS curve corresponds to the step of the SM which follows the state with E = Es, as de- 
scribed above. 

^ Of these two intervals, the relevant one is that in which the derivative dE/dt in the nonlinear 
resonance equations (see Eq. (16) below) is positive or negative, for the case of the upper or lower 
GSS curve respectively. 

' I.e. E determined by Eq. (7) for any \j/ except from the vicinity of multiples of 7t. As shown in 
1431 . Eq. (7) is irrelevant to the boundary of the chaotic layer in the range of iff close to multiples 
of 7t while the boundary in this range of \jf still lies in the resonance range of energies, where 
m{E) « £o(±). 
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slow angle ^ slow angle tjjr slow angle 



Fig. 1 Schematic figure illustrating the formation of the peak of the function A £"„„ {(Of): (a) Wf = 
fflmax; (b) (Of < cOmax; (c) (Of > OJmax- The relevant (lower) GSS curve is shown by the dotted 
line. The relevant trajectories of the resonance Hamiltonian are shown by solid lines. The lower 
boundary of the layer is marked by a thick solid line: in (a) and (b) the lower boundary is formed 
by the lower part of the resonance separatrix while, in (c) it is formed by the resonance trajectory 
tangent to the GSS curve. The dashed line mai'ks, for a given (Of, the maximal deviation of the 
lower boundary from the separatrix energy E^. 



self-intersection) touches or intersects the upper/lower GSS curve while the up- 
per/lower part does not, then the upper/lower boundary of the layer is formed by 
the upper/lower part of this self-intersecting trajectory (Figs. 1(a) and 1(b)); (ii) oth- 
erwise the boundary is formed by the resonance trajectory tangent to the GSS curve 
(Fig. 1(c)). It is shown below that, in both cases, the variation of the energy along 
the resonance trajectory is larger than the separatrix split 5 by a logarithmically 
large factor o= ln(l//!). Therefore, over the boundary of the SM chaotic layer the 

largest deviation of the energy from the separatrix value, AEi^\ may be taken, in 
the leading-order approximation, to be equal to the largest variation of the energy 
along the resonance trajectory forming the boundary, while the latter trajectory can 
be entirely described within the resonance Hamiltonian formalism. 

Finally, we mention that, as is obvious from the above description of the bound- 
ary, ziiiim ' = AEim\(Of) possesses a local maximum AE^l^sm at (Of for which the 
resonance separatrix just touches the corresponding GSS curve (see Fig. 1(a)). 



3 Single-Separatrix Chaotic Layer 

It is clear from Sec. 2 above that AE^x.sm is equal in leading order to the width 
AEnr of the nonlinear resonance which touches the separatrix. In Sec. 3.1 below, 
we roughly estimate AEf^R in order to classify two different types of systems. Sees. 
3.2 and 3.3 present the accurate leading-order asymptotic theory for the two types 
of systems. The next-order correction is estimated in Sec. 3.4, while a discussion is 
presented in Sec. 3.5. 
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3.1 Rough estimates. Classification of systems. 

Let us roughly estimate AE^r: it will turn out that it is thus possible to classify all 
systems into two different types. With this aim, we expand the perturbation V into 
two Fourier series in t and in i// respectively: 

y = ^£y(''(£,V/)e-""^'+cx.= i^yf (£)e'(^-'^-''«/')+c.c. (8) 

^ / ^ l,k 

As in standard nonlinear resonance theory lfT0ll23l 155115111521 . we single out the 
relevant (for a given peak) numbers K and L for the blind indices k and / respectively, 
and denote the absolute value of V^'^' as Vq: 

Vo{E)^\vII'\e)\. (9) 

To estimate the width of the resonance roughly, we use the pendulum approxi- 
mation of resonance dynamics lITOl l23l l55l ISTl l52l [Tl : 



AEnr - ^%hVQ(Of/\A(o/AE\. (10) 

This approximation assumes constancy of d(o/dE in the resonance range of ener- 
gies, which is not the case here; in reality, (o{E)°^ \/\vi{\/\E — Es\) in the vicinity of 
the separatrix l53[l0l|23|55l|5ll|52l[I]|2l[3l|43[35lB^ that the relevant 

derivative Idco/disl ~ {(Or^'')'^ / {(Oq\E — Es\) varies strongly within the resonance 
range. However, for our rough estimate we may use the maximal value of \E —Es\, 
which is approximately equal to AEnr. If (Of is of the order of o,-*' ^ a)t)/ln(l /h), 
then Eq. (10) reduces to the following approximate asymptotic equation for AEnr'- 

AEnr^Vo{E =E,±AENR)li\n{l/h), (11) 

The asymptotic solution of Eq. (11) depends on Vo(£'.s ± AEnr) as a function of 
AEnr- In this context, all systems can be divided in two types. 

I The separatrix of the unperturbed system has two or more saddles while the rele- 
vant Fourier coefficient y'^' = V {E, \ff) possesses different values on adjacent 
saddles. Given that, for E Eg, the system stays most of time near one of the 
saddles, the coefficient — > Es,\l/) as a function of \j/ is nearly a "square 
wave": it oscillates between the values at the different saddles. The relevant K 
is typically odd and, therefore, Vo{E Es) approaches a well defined non-zero 
value. Thus, the quantity Vo{E = Es ± AEnr) in Eq. (1 1) may be approximated 
by this non-zero limit, and we conclude therefore that 

AEnr°- h\n{l/h), h^O. (12) 

n Either (i) the separatrix of the unperturbed system has a single saddle, or (ii) it 
has more than one saddle but the perturbation coefficient V'^' is identical for all 
saddles. Then y(^) (£ i//), as a periodic function of y/, significantly differs 
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from its value at the saddle(s) only during a small part of the period in y/: this 
part is ~ co(£)/coo 1 / ln(l /|£, I ). Hence, yo(£, ±zi£A,fl) 1/ ln(l /zi£A,x). 
Substituting this value in Eq. (1 1), we conclude that 

AEm^h, h^O. (13) 

Thus, for systems of type I, the maximum width of the SM chaotic layer is propor- 
tional to li times a logarithmically large factor 0= ln(l /li) while, for systems of type 
II, it is proportional to h times a numerical factor 

As shown below, the variation of energy in between the instants relevant to the 
SM is ~ h, i.e. much less than AEnr (12) for systems of the type I, and of the same 
order as AEnr (13) for systems of type II. Therefore, one may expect that the max- 
imum width of the layer for the original Hamiltonian system (1), AE^^\ is at least 
roughly approximated by that for the SM, zi^v,^', so that the above classification of 
systems is relevant to ZiZs'^*' too. This is confirmed both by numerical integration 
of the equations of motion of the original Hamiltonian system and by the accurate 
theory presented in the next two sub-sections. 



3.2 Asymptotic theory for systems of type I. 

For the sake of clarity, we consider a particular example of a type I system; its 
generalization is straightforward. 

We choose an archetypal example: the ac-driven pendulum (sometimes referred 
to as a pendulum subject to a dipole time-periodic perturbation) 1551 l42l |45]| : 



H = HQ + hV, (14) 
//o = — -cos(^), V = -qcos{cOft), /!<1. 

Fig. 2 presents the results of numerical simulations for a few values of h and sev- 
eral values of o)/. It shows that: (i) that the function AE^^\a)f) indeed possesses 
sharp peaks whose heights greatly exceed the estimates by the heuristic 1551 . adia- 
batic IT4I and moderate-frequency l29l theories (see inset); (ii) as predicted by our 
rough estimates of Sec. 3.1, the 1st peak of AE^^^^cOf) shifts to smaller values of 
(Of while its magnitude grows, as h decreases. Below, we develop a leading-order 
asymptotic theory, in which the parameter of smallness is l/ln(l /h), and compare 
it with results of the simulations. 

Before moving on, we note that the SM (approximated in the relevant case by 
nonlinear resonance dynamics) considers states of the system only at discrete in- 
stants. Apart from the variation of energy within the SM dynamics, a variation of 
energy in the Hamiltonian system also occurs in between the instants relevant to the 
SM. Given that O/ <C 1, this latter variation may be considered in adiabatic approxi- 
mation and it is of the order of h lT4l[35]| . It follows from the above rough estimates. 
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and from the more accurate consideration below, that the variation of energy within 
the SM dynamics for systems of type I is logarithmically larger i.e. larger by the 
factor ln(l /h). The variation of energy in between the instants relevant to the SM 
may therefore be neglected to leading-order for systems of type I: AE^^'^ ~ AEsm^ . 
For the sake of notational compactness, we shall henceforth omit the subscript "sm" 
in this subsection. 

For the system (14), the separatrix energy is equal to 1, while the asymptotic (for 
E Es) dependence (£i{E) is 1551 : 

"^^)- ln(32/|I.-^|) ' '''' 
£, = 1, |£,-£|<1. 

Let us consider the range of energies below E^ (the range above E^ may be con- 
sidered in an analogous way) and assume that (Of is close to an odd multiple of 

(Or \ The nonlinear resonance dynamics of the slow variables in the range of 




Fig. 2 Computer simulations for the ac driven pendulum (14) (an archetypal example of type 1): 
the deviation of the lower boundary of the chaotic layer from the separatrix, normalized by 

the perturbation amplitude h, is plotted as a function of the perturbation frequency (Of, for vari- 
ous h. The inset presents the same data but with a logarithmic ordinate and with the estimates by 
the heuristic 1551 . adiabatic 1141 and moderate-frequency 1291 theories. The heuristic estimate is 
shown by the dotted line: as an example of the heuristic estimate, we use the formula from 1551 : 
AE^-^^h = iTKOf / co&\i{K(£if /2) . The adiabatic and moderate-frequency estimates are shown by 
the dashed line: the adiabatic estimate for 4£'^'(ft)/) is equal approximately to 2k; the estimate 
following from the results of the work (29| for (Of ~ 1 is of the same order, so that it is schemati- 
cally represented in the inset in Fig. 2 by the same line as for the adiabatic estimate (dashed line). 
The inset shows explicitly that the simulation results exceed the estimates of the former theories 
by 1 or 2 orders of magnitude, over a wide range of frequencies. 
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approximately resonant energies may be described as follows ||43l 1411 (cf. also 
HSIlllllllSlllllII]): 



dt df ' df dl ' ^ ^ 

H{I,\ir)= [ dl {nco- COf) - nhq„cos{\if) 

Jl(Es) 

= n{E - Es) — (Of{I — I{Es)) - nhq„cos{\jf) , 

i-E ^ 

I^I{E)= — ^, E = Ho{p,q), 

= COft, 

fl dq 
W^7t + signip)coiE) \ .+27tl, 

JimJE) ^2[E-U{q)) 

1 /■2'i 

q„ = q„{E)^— dv/^(£',r)cos(«V/), 
271 Jq 

\nCO- (Of\ CO, n = 2j-l, j= 1,2,3,..., 

where / and y/ are the canonical variables action and angle respectively ifTOl |23l l55l 
|5Tl|52][l]; iiniin is the minimal energy over all q,p, E = Ho{p,q); qmm{E) is the 
minimum coordinate of the conservative motion with a given value of energy E; I is 
the number of right turning points in the trajectory [^(t)] of the conservative motion 
with energy E and given initial state (qo^po). 

The resonance Hamiltonian H{I, y) is obtained in the following way. First, the 
original Hamiltonian H is transformed to action-angle variables I — y/. Then it is 
multiplied by « and the term (Of I is extracted (the latter two operations correspond 
to the transformation yf ^ \if = n\j/~ CO ft). Finally, the result is being averaged over 
time i.e. only the resonance term in the double Fourier expansion of the perturbation 
is kept (it may be done since the effect of the fast-oscillating terms on the dynamics 
of slow variables is small: see the estimate of the corrections in Sec. 3.4 below). 

Let us derive asymptotic expression for /(£), substituting the asymptotic expres- 
sion (15) for Co{E) into the definition of I{E) (16) and carrying out the integration: 

'm-'(«-^(>n(5^) + l). (17) 

As for the asymptotic value q„{E E^), it can be seen that q{E ^ E^, Xj/), as a 
function of y/, asymptotically approaches a "square wave", oscillating between —n 
and TT, so that, for sufficiently small j, 



qy_,{E^Es)c^ {-!}'+' (18) 

=0, 
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The next issue is the analysis of the phase space of the resonant Hamiltonian 
(16). Substituting ft (16) into the equations of motion (16), it can be seen that their 
stationary points have the following values of the slow angle 

V/+ = 7r, v/_=0, (19) 
while the corresponding action is determined by the equation 

Aq„ 

n(0- COf^nh— =Q, n = 2y-l, (20) 

where the sign "=F"corresponds to (19). 

The term °^ h in (20) may be neglected to leading-order (cf. ifTOl l23l l55l ISTl l52l 
[Tl l43ll4ll ). and Eq. (20) reduces to the resonance condition 

{2j-\)(0{El-'^) = (0f, (21) 
the lowest-order solution of which is 

£.-£(^')^32exp(-fcil^). (22) 

Eqs. (19) and (22) together with (17) explicitly determine the elliptic and hyper- 
bolic points of the Hamiltonian (16). The hyperbolic point is often referred to as a 
"saddle" and corresponds to or in (19) for even or odd j respectively. The 
saddle point generates the resonance separatrix. Using the asymptotic relations (17) 
and (18), we find that the resonance Hamiltonian (16) takes the following asymp- 
totic value in the saddle: 



Hsaddle — % — 2/l 

^'^32.J-'^^^\-2h. (23) 
n \ cof J 

The second asymptotic equality in (23) takes into account the relation (22). 

As explained in Sec. 2 above, AE'^^^cOf) possesses a local maximum at cOf for 
which the resonance separatrix is tangent to the lower GSS curve (Fig. 1(a)). For the 
relevant frequency range O/ 0, the separatrix split (which represents the maxi- 
mum deviation of the energy along the GSS curve from Es) approaches the following 
value ll55l in the asymptotic limit h-^Q 

8~2nh, a)/<l. (24) 

As shown below, the variation of energy along the relevant resonance trajectories 
is much larger Therefore, in the leading-order approximation, the GSS curve may 
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simply be replaced by the separatrix of the unperturbed system i.e. by the horizontal 
line E = Eg or, equivalently, / = I (Eg). Then the tangency occurs at iff, shifted from 
the saddle by n, so that the condition of tangency is written as 

H,acUlle =H{I= /(£,), Y = ^saddle + T^) ^ 2h. (25) 

Substituting here Hsaddie (23), we finally obtain the following transcendental 
equation for COmax- 

„xpW.?efIi, ..2/^. (26, 

Fig. 3(b) demonstrates the excellent agreement between Eq. (26) and simulations of 
the Hamiltonian system over a wide range of h. 

In the asymptotic limit /; 0, the lowest-order explicit solution of Eq. (26) is 





Fig. 3 An archetypal example of a type 1 system: the ac-driven pendulum (14). Comparison of the- 
ory (solid lines) and siinulations (circles) for: (a) the deviation AE^^'>(af) of the lower boundary 
of the chaotic layer from the separatrix, normalized by the perturbation amplitude h, as a function 
of the perturbation frequency (Of, for h = 10"''; the theory is from Eqs. (26), (31), (32), (38), (39) 
and (41) (note the discontinuous drop by the factor e from the maximum to the right wing), (b) 
The frequency of the 1st maximum in AE^^mf) as a function of h; the theory is from Eq. (26). 
(c) The 1st maximum in 4£("' {mf) jh as a function of h; the theory is from Eqs. (34) and (26). 
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C^kaU jA^^L_^, ;.l,2,...«ln(-j. (27) 

As follows from Eq. (26), the value of — Er^'^ (22) for cOf = fi),14x is 

E,-Ey\(Of = o^£) = ^. (28) 
Its leading-order expression is: 

E. Ei.^\cOf = co^^ ^ ^ in (^^) , - 0. (29) 

If (Of < CoHlx then, in the chaotic layer, the largest deviation of energy from the 

separatrix value corresponds to the minimum energy E^l^ on the nonlinear reso- 
nance separatrix (Fig. l(a,b)), which occurs at xjr shifted by n from the saddle. The 
condition of equality of H at the saddle and at the minimum of the resonance sepa- 
ratrix is written as 

H,addle = HiliE^^l), f^addle + 7t). (30) 

Let us seek its asymptotic solution in the form 



E,-E± = AeI'^ = il+y)iE,-El''^)c^il+y)32exp( ""^^^ ^' 



(Of 

y^l. (31) 

Substituting (31) and (23) into Eq. (30), we obtain fory the following transcen- 
dental equation: 



(l+y)ln(l+y)-y=^^^j^x/exp(x/), (32) 
Xf=^^ (Of<(O^L, y>0, 

(Of 

where onl^ is given by Eq. (26). 

Eqs. (31) and (32) describe the left wing of the j-th peak of AE^-^"^ (®f )■ Fig- 3(a) 
demonstrates the good agreement between our analytic theory and simulations for 
the Hamiltonian system. 

It follows from Eq. (26) that Eq. (32) for cOf = CO^lx reduces to the relation 
ln(l +y) = 1, i.e. 

l+y(cOiL4l)=e. (33) 
It follows from Eqs. (33), (31) and (28) that the maximum for a given peak is: 
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AeIII^^E,~E±{co^±)^^-^. (34) 

Fig. 3(c) shows the excellent agreement of this expression with our simulations 
of the Hamiltonian system over a wide range of h. 
The leading-order expression for is: 

~ _l£^Lin(8(2y - l)//i), ft ^ 0, (35) 

which confirms the rough estimate (12). 

As (Of decreases, it follows from Eq. (32) that y increases exponentially sharply. 

In order to understand how Ae\'^'' decreases upon decreasing (Of, it is convenient to 
rewrite Eq. (31) re-expressing the exponent by means of Eq. (32): 

^^""■'/' = ».,|ln(l+.)-,/(l+.)) - 

It follows from Eqs. (32) and (36) that Ae\^^ decreases power-law-like when (Of is 

decreased. In particular, Ae\^^ °^ \/ {(0^1^ — (Of) at the far part of the wing. 

As for the right wing of the peak, i.e. for (Of > OJmax, over the chaotic layer, the 
largest deviation of energy from the separatrix value corresponds to the minimum 
of the resonance trajectory tangent to the GSS curve (Fig. 1(c)). The value of iff at 
the minimum coincides with "^saddle- In the leading-order approximation, the GSS 
curve may be replaced by the horizontal line / = I{Es), so that the tangency occurs 

at xjf = ^fsaddie + ^- Then the energy at the minimum E^J^^^ can be found from the 
equation 

//(/(£,), V^W.//. + 7t)^ HiI{E2l), Saddle) (37) 

Let us seek its asymptotic solution in the form 



E, - ^ AEj:^^ = z{E, - c z32 exp ^^^^^Jl 
0<z<l, (38) 
Substituting (38) into (37), we obtain for z the following transcendental equation: 



z(l+ln(l/z)) = ^^^^x/exp(x/) (39) 
xf=^'^ — 11^ to/>a)i4x, 0<z<l, 

(Of 
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where Oinax is given by Eq. (26). Eqs. (38) and (39) describe the right wing of the 
7-th peak of AE''^\a)f). Fig. 3(a) demonstrates the good agreement between our 
analytic theory and simulations. 

It follows from Eq. (26) that the solution of Eq. (39) for cOf Co!nlx is z ^ 1, 
so the right wing starts from the value given by Eq. (28) (or, approximately, by Eq. 
(29)). Expressing the exponent in (38) from (39), we obtain the following equation 

a)/(l+ln(l/z)) 

It follows from Eqs. (39) and (40) that ZiiSr^' decreases power-law-like for increas- 
ing (Of. In particular, ZiZ?,-^' 1/ (o)/ — CO^lx) in the far part of the wing. Further 
analysis of the asymptotic shape of the peak is presented in Sec. 3.5 below. 

Beyond the peaks, the function AE^^\a)f) is logarithmically small in compari- 
son with the maxima of the peaks. The functions AE^''\cOf) and AEr^\(i)f) in the 
ranges beyond the peaks are also logarithmically small. Hence, nearly any function 
of AEr'\cOf) and AEj-''^^\cOf) which is close to AEi^\(Of) in the vicinity of (O^^x 
and to AEj''^^\a)f) in the vicinity of fi)iS4x'' while being sufficiently small beyond 
the peaks may be considered as an approximation of the function AE^^'f (o/) to log- 
arithmic accuracy with respect to the maxima of the peaks, 4-Emax and ZiEmalt'^ in 
the whole range [co^L,(okiL^^]- One of the easiest options is the following: 



AE^-\(Of) = AEj^\(Of) for (Of < coial, 

AE'^-\(Of) = max{AEj^^\cOf),AEl^+^\cOf)} for cojH < cOf < coi^i+'' 



ax 5 



7 = l,2,...«^j^. (41) 

ax 

We used this function in Fig. 3(a), and the analogous one will also be used in the 
other cases. 

In fact, the theory may be generahzed in such a way that Eq. (41) would approxi- 
mate AE^'^^cOf) well in the ranges far beyond the peaks with logarithmic accuracy, 

even with respect to AE^^\(Of) itself rather than to AEmlx only (cf. the next sec- 
tion). However, we do not do this in the present case, being interested primarily in 
the leading-order description of the peaks. 

Finally, we demonstrate in Fig. 4 that the lowest-order theory describes the 
boundary of the layers quite well, even in the Poincare section rather than only 
in energy/action. 
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Fig. 4 Some characteristic Poincare sections in the 2;r-interval of the energy-angle plane for the 
system (14) with h = lO"' and (Of equal to: (a) 0.236 (maximum), (b) 0.21 (left wing), (c) 0.25 
(right wing). Results of the numerical integration of the equations of motion for the original Hamil- 
tonian (14) are shown by (red) dots. The NR separatrix calculated in the leading-order approxima- 
tion (i.e. by the integration of the resonant equations of motion (16) in which a>{E), 1(E) and 
q\ (E) are approximated by the explicit formulae (15), (17) and (18) respectively) is drawn by the 
(black) solid line. The NR trajectory (calculated in the leading-order approximation) tangent to the 
line E = £, is drawn by the (blue) dashed line. The outer boundary (marked by a thicker line) is 
approximated by: the lower part of the NR separatrix in cases (a) and (b), and by the tangent NR 
trajectory in case (c) The boundary of the island of stability in the cases (a) and (b) is approximated 
by the tangent NR trajectory (which coincides in the case (a) with the NR separatrix). 



3.3 Asymptotic theory for systems of type II. 

We consider two characteristic examples of type II systems, corresponding to the 
classification given in Sec. 3.1. As an example of a system where the separatrix of 
the unperturbed system possesses a single saddle, we consider an ac-driven Duffing 
oscillator |[Tl[T8]|29]|40l. As an example of the system where the separatrix possesses 
more than one saddle, while the perturbation takes equal values at the saddles, we 
consider a pendulum with an oscillating suspension point lU [18] |29l [34l [35l . The 
treatment of these cases is similar in many respects to that presented in Sec. 3.2 
above. So, we present it in less detail, emphasizing the differences. 
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3.3.1 AC-driven Duffing oscillator. 

Consider the following archetypal Hamiltonian 

H = Hq + IiV, (42) 

2 2 4 

//o = y-y + ^, V = -qcos{cOft), /i< 1. 

The asymptotic dependence of Co{E) on E for E below the separatrix energy 
£v = is the following ifTlfTSl 



£, = 0, 1- 

Correspondingly, the resonance values of energies (determined by the condition 
analogous to (21)) are 

E,-EjJ^ = 16exp (-^] , j = 1,2,3, ... (44) 



COf 

The asymptotic dependence of I{E) is 

/W-/K.)-^(ln(^) + l). (45, 

The nonlinear resonance dynamics is described by the resonance Hamiltonian H 
which is identical in form to Eq. (16). Obviously, the actual dependences co{E) and 
I{E) are given by Eq. (43) and (45) respectively. The most important difference is 
in qj{E): instead of a non-zero value (see (18)), it approaches as £ ^ E^. Namely, 
itis°c to(£) lUlini: 

<!jiE)-^co{E), ; = l,2,...«^, (46) 

i.e. qj is much smaller than in systems of type I (cf. (18)). Due to this, the resonance 
is "weaker". At the same time, the separatrix split 5 is also smaller, namely ^ licOf 
(cf. Il43l ) rather than ~ li as for the systems of type 1. That is why the separatrix 
chaotic layer is still dominated by resonance dynamics while the matching of the 
separatrix map and nonlinear resonance dynamics is still valid in the asymptotic 
limit /j^OlEl ■ 

Similarly to the previous section, we find the value of H in the saddle in the 
leading-order approximatior(§ 



* The only essential difference is that q„ at the saddle is described by Eq. (46) rather than by Eq. 
(18). 
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where Es — is given in (44). 

As before, the maximum width of the layer corresponds to (Of, for which the 
resonance separatrix is tangent to the GSS curve (Fig. 1(a)). It can be shown B3l that 
the angle of tangency asymptotically approaches iffsaddie + = tt while the energy 
still lies in the resonance range. Here Co{E) w ©r "^^(Of/ j. Using the expressions 
for H{E,\jf) (cf. (16)), I{E) (45), qj{E) (46), and taking into account that in the 
tangency E < 8 ^ liCOf <C h, to leading-order the value of H at the tangency reads 

h 

Htangency — 0)f--j=. (48) 

Allowing for Eqs. (47) and (48), the condition for the maximum, Hsuddie — 

Htangency, reduces tO 

E,-Ei'\a)^l)c^27tV2h. (49) 

Thus these values E^ — /ir^' are logarithmically smaller than the corresponding 
values (28) for systems of type 1. 

The values of (Of corresponding to the maxima of the peaks in AE'^^^cOf) are 
readily obtained from (49) and (44): 

The derivation to leading order of the shape of the peaks for the chaotic layer 
of the separatrix map, i.e. within the nonlinear resonance (NR) approximation, is 
similar to that for type 1. So, we present only the results, marking them with the 
subscript "NR". 

The left wing of the Jth peak of AE''^J{(Of) is described by the function 

- .6(1 +,)exp , ,51, 



(Of J ln(l+y)-3;/(l+3;)' 



(Of < co^L 



where y is the positive solution of the transcendental equation 

7th / In j 



(l+y)ln(l+y)-3'=^exp(^-^j, y>0. (52) 
In common with the type I case, 1 +y{(0^]ix) = so that 
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^EliLm = - Ej:'\c4il)) ~ 2;re^/2/1. (53) 

Eq. (53) confirms the rough estimate (13). The right wing of the peak is described 
by the function 



^4i(a.,)^lfeexpf-^).^i^^, (54) 



(Of J l+ln(l/z)' 



(Of > CO^lx, 



where z < 1 is the solution of the transcendental equation 

z(l+ln(l/z)) = ^exp(^^), 0<z<l. (55) 

As in the type I case, z{cOf — > fi)m2t) 1. 

It follows from Eqs. (49) and (53) that the typical variation of energy within 
the nonlinear resonance dynamics (that approximates the separatrix map dynamics) 
is cx h. For the Hamiltonian system, the variation of energy in between the dis- 
crete instants corresponding to the separatrix map 1551 ISTl l52l [l] |43] [3T| is also 
h. Therefore, unlike the type 1 case, one needs to take it into account even at the 
leading-order approximation. Let us consider the right well of the Duffing potential 
(the results for the left well are identical), and denote by f^^ the instant at which the 
energy £ at a given k-th step of the separatrix map is taken: it corresponds to the 
beginning of the A:-th pulse of velocity ||55] |43l i.e. the corresponding q is close to 
a left turning point qifp in the trajectory [^'(t)]. Let us also take into account that 
the relevant frequencies are small so that the adiabatic approximation may be used. 
Thus, the change of energy from fj^ up to a given instant f during the following pulse 
of velocity {t — tj^^ 1) may be calculated as 



AE ^ / dTqhcos{a)fT) hcos{cOftii) / dxq 

= hcos{(Oftk){q{t)-qitp) (56) 

For the motion near the separatrix, the velocity pulse corresponds approximately 
to i/A = (see the definition of (16)). Thus, the corresponding slow angle is xfr = 

i\l/-(Oftk~-(Oftk. 

For the left wing of the peak of AE'^^^cOf) (including also the maximum of 
the peak), the boundary of the chaotic layer of the separatrix map is formed by the 
lower part of the NR separatrix. The minimum energy along this separatrix occurs at 
iff = 7t. Taking this into account, and also noting that \fr ~ —cOftk, we conclude that 
cos{cOftii) ~ — 1. So, AE < 0, i.e. it does lower the minimum energy of the layer of 
the Hamiltonian system. The maximum reduction occurs at the right turning point 

9rtp- 
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max{\AE\) ~ h{qr,p - qi,p) = Vlh. (57) 
We conclude that the left wing of the y-th peak is described as follows: 

AE\'\(Of)c^AE\'^^{(Of) + ^h, cOf<co^L (58) 

where AE\^^^{cOf) is given by Eqs. (51)-(52). In particular, the maximum of the 
peak is: 

AeHI, ~ (2;re + l)^/2h w 25.6/?. (59) 

For the right wing of the peak, the minimum energy of the layer of the separatrix 
map occurs when ^r coincides with saddle (Fig- 1(c)) i.e. is equal to 0. As a result, 
co?,{(i)ftii) ~ 1 and, hence, AE > 0. So, this variation cannot lower the minimum 

energy of the layer for the main part of the wing, i.e. for cOf < coj/],,^! where ColH^j is 

defined by the condition AE^.^^j^ = max{\AE\) = ^/ih. For (Of > (olHj, the minimal 
energy in the layer occurs at\j/=n, and it is determined exclusively by the variation 
of energy during the velocity pulse (the NR contribution is close to zero at such iff). 
Thus, we conclude that there is a bending of the wing at cOf = O/^L/' 



AEi^\cOf) = AE%{cOf), o^2x <C0f< coild^ 



AE9\cOf) = V2h, Of > (ol^'^ 



ln(8V2//7) + l - 2;r' 



bend ' 

(60) 



where AE^J^i^{(Of) is given by Eqs. (54) and (55). 

Analogously to the previous case, AE'^^^aj) may be approximated over the 

whole frequency range by Eq. (41) with Ae\^^ and Ae'"/^ given by Eqs. (58) and 
(60) respectively. Moreover, unlike the previous case, the theory also describes ac- 
curately the range far beyond the peaks: ZiE^^' is dominated in this range by the 
velocity pulse contribution AE, which is accurately taken into account both by Eqs. 
(58) and (60). 

Fig. 5 shows very reasonable agreement between the theory and simulations, 
especially for the 1st pealfl 



' The disagreement between theory and simulations for the magnitude of the 2nd peak is about 
three times larger than that for the 1st peak, so that the height of the 2nd peak is about 30% smaller 
than that calculated from the asymptotic theory. This occurs because, for the energies relevant to 
the 2nd peak, the deviation from the separatrix is much higher than that for the 1st peak. Due to 
the latter, the Fourier coefficient q2{E) for the relevant E is significantly smaller than that obtained 
from the asyinptotic formula (42). In addition, the velocity pulse contribution AE also significantly 
decreases while the separatrix split increases as (Of becomes ~ 1. 
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3.3.2 Pendulum with an oscillating suspension point 

Consider the archetypal Hamiltonian ifTI [TSl |29l l34l |351 



2 

//() = ^ +cos(^), V = ~cos,{q)cos,{(Oft), (61) 

Though the treatment is similar to that used in the previous case, there are also 
characteristic differences. One of them is the following: although the resonance 
Hamiltonian is similar to the Hamiltonian (16), instead of the Fourier component 
of the coordinate, q„, there should be the Fourier component of cos,{q), V,,, which 
can be shown to be: 




Fig. 5 An archetypal example of a type II system: the ac driven Duffing oscillator (42). Com- 
parison of theory (solid lines) and simulations (circles): (a) the deviation AE^^\(Of) of the lower 
boundary of the chaotic layer from the separatrix, normalized by the perturbation amplitude h, as 
a function of the perturbation frequency (Of, for /; = 10^^; the theory is from Eqs. (41), (50), (51), 
(52), (54), (55), (58) and (60) (note the discontinuous drop from the maximum to the right wing); 

(b) the frequency of the 1st maximum in (co^ ) as a function of h; the theory is from Eq. (50); 

(c) the 1 St maximum in (tt"/ )/ h as a function of /;; the theory is from Eq. (59). 
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V2jC^i-iy+'^coiE), (62) 

^2,-1=0, 

In 1 

2n Jo 

The description of the chaotic layer of the separatrix map at the lowest order, i.e. 
within the NR approximation, is similar to that for the ac-driven Duffing oscillator 
So we present only the results, marking them with the subscript ''NR". 

The frequency of the maximum of a given j-th peak is: 

cokL - ^^^y j = 1,2, ... « ln(4//0. (63) 

This expression agrees well with simulations for the Hamiltonian system (Fig. 6(b)). 
To logarithmic accuracy, Eq. (63) coincides with the formula following from Eq. 
(8) of ||34l (reproduced in ||35l as Eq. (21)) taken in the asymptotic limit li ^ 
(or, equivalently, a),nax ^ 0). However, the numerical factor in the argument of the 



logarithm in the asymptotic formula following from the result of Il34[|35l is half our 
value: this is because the nonlinear resonance is approximated in ll34l [35l by the 
conventional pendulum model which is not valid near the separatrix (cf. our Sec. 
3.1 above). 

The left wing of the jth peak of Asj^j^ {(Of) is described by the function 

44^(0),) = 32(1 +y)exp f-^) ^ , , (64) 



(Of) ln(l+3;)-y/(l+3;)- 



(Of < (okL, 



where y is the positive solution of the transcendental equation 

(l+3')ln(l+3')-3'= Jexp(^^^, y>0. (65) 
Similarly to the previous cases, 1 +y{(Omlx) ~ e. Hence, 

^eI^Lnr = - (o^i)) = 8e/j. (66) 

Eq. (66) confirms the rough estimate (13). The right wing of the peak is described 
by the function 



<.K'„ UnA = ^?7pxn = 



AE-,(co,) = 32zexp ^-^ j . (67) 

(Of > (okL, 

where z < 1 is the solution of the transcendental equation 
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z(l+ln(l/z)) = %xp( 0<z<l. (68) 

4 \0} f J 

Similarly to the previous cases, z{cOf CO^L) 1- 

Now consider the variation of energy during a velocity pulse. Though the final 
result looks quite similar to the case with a single saddle, its derivation has some 
characteristic differences, and we present it in detail. Unlike the case with a single 
saddle, the pulse may start close to either the left or the right turning point, and the 
sign of the velocity in such pulses is opposite 1551 |43l . The angle y/ in the pulse is 
close to —n/2 or n/2 respectively. So, let us calculate the change of energy from 
the beginning of the pulse, f^, until a given instant t within the pulse: 

AE = — f drqlidV /dq = h f dTq{—sm{q)cos{cOfT)) 

~/icos(a)/fA.) / dT^(-sin(^)) ~ /zcos(a)/fA.)(cos(^(f)) - 1). (69) 

Here, the third equality assumes adiabaticity while the last equality takes into ac- 
count that the turning points are close to the maxima of the potential i.e. close to a 
multiple of 2n (where the cosine is equal to 1). 

The quantity AE (69) takes its maximal absolute value at ^ = TT. So, we shall 
further consider 



^^max = -2/icos(a)/fA.) = -2/icos(2;>i - ft) = (-l)^'+'2/icos(v>'0. (70) 

The last equality takes into account that, as mentioned above, the relevant y/^. is 
either — 7r/2 or k/2. For the left wing, the value of f at which the chaotic layer of 
the separatrix map possesses a minimal energy corresponds to the minimum of the 
resonance separatrix. It is equal to TT or if the Fourier coefficient V2j is positive or 
negative, i.e. for odd or even j, respectively: see Eq. (63). Thus AEma\ = —2h for 
any j and, therefore, it does lower the minimal energy of the boundary. We conclude 
that 

AE\'\(Of) ~ AE\'li,{(Of) + 2h, (Of < a^L (71) 

where AEj''^j^{cOf) is given by Eqs. (64)-(65). In particular, the maximum of the 
peak is: 

^E^L =i (4e + 1 )2h w 23.7/1. (72) 

The expression (72) confirms the rough estimate (13) and agrees well with sim- 
ulations (Fig. 6(c)). At the same time, it differs from the formula which can be ob- 
tained from Eq. (10) of El (using also Eqs. (1), (3), (8), (9) of 01) in the asymp- 
totic Umit /j — > 0: the latter gives for AEj-il,^ the asymptotic value 32/j. Though this 
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result 11341 (referred to also in ||35l ) provides for the correct functional dependence 
on h, it is quantitatively incorrect because (i) it is based on the pendulum approxi- 
mation of the nonlinear resonance while this approximation is invalid in the vicinity 
of the separatrix (see the discussion of this issue in Sec. 3.1 above), and (ii) it does 
not take into account the variation of energy during the velocity pulse. 

The right wing, by analogy to the case of the Duffing oscillator, possesses a 

bend at cOf ~ (Ohemj where ^/i^jij^ = jziiiniaxl = 2/z, corresponding to the shift of the 
relevant iff for n. We conclude that: 




10 - 



10"' 10"' 10"'' 10""' 10"" 

h 

Fig. 6 An archetypal example of a type II system: the penduluin with an oscillating suspen- 
sion point (61). Comparison of theory (solid lines) and simulations (circles): (a) The deviation 
AE^^''{o)f) of the lower boundary of the chaotic layer from the separatrix, normalized by the per- 
turbation amplitude h, as a function of the perturbation frequency (Df, for h = 10^^; the theory is 
by Eqs. (41), (63), (64), (65), (67), (68), (71) and (73) (note the dicontinuous drop from the max- 
imum to the right wing), (b) The frequency of the 1st maximum in AE^^O)f) as a function of 
the theory is from Eq. (63). (c) The 1st maximum in ziS'^' (ft)/) /h as a function of h; the theory is 
from Eq. (72). 
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CO. 



U) _ 



2%j 



(73) 



'bend 



ln(16//7)-3' 



where ^^^^^(tO/) is given by Eqs. (66) and (67). 

Similarly to the previous case, both the peaks and the frequency ranges far be- 
yond the peaks are well approximated by Eq. (41), with AE^^^^ and Zi/ir^^ given by 
Eqs. (71) and (73) respectively (Fig. 6(a)). 

3.4 Estimate of the next-order corrections 

We have calculated explicitly only the leading term AE in the asymptotic expansion 
of the chaotic layer width. Explicit calculation of the next-order term AE^"'^'^^ is 
possible, but it is rather complicated and cumbersome: cf. the closely related case 
with two separatrices 1431 (see also Sec. 4 below). In the present section, where the 
perturbation amplitude h in the numerical examples is 4 orders of magnitude smaller 
than that in B3l . there is no particular need to calculate the next-order correction 
quantitatively. Let us estimate it, however, in order to demonstrate that its ratio to 
the leading term does vanish in the asymptotic limit h 0. 

( fje Xt ) 

We shall consider separately the contribution 4 £■„■ stemming from the various 
corrections within the resonance approximation (16) and the contribution AE^""^^^ 
stemming from the corrections to the resonance approximation. 

The former contribution may be estimated in a similar way to the case considered 
in ||43l : it stems, in particular, from the deviation of the GSS curve from the separa- 
trix (this deviation reaches 8 at certain angles: see Eq. (7)) and from the difference 
between the exact resonance condition (20) and the approximate one (21). It can be 
shown that the absolute value of the ratio between zif,!""'* and the leading term is 
logarithmically small (cf. 1431 ): 



Let us turn to the analysis of the contribution AE, , i.e. the contribution stem- 
ming from the corrections to the resonance Hamiltonian (16). It is convenient to 
consider separately the cases of the left and right wings of the peak. 

As described in Sees. 3.2 and 3.3 above, the left wing corresponds in the leading- 
order approximation to formation of the boundary of the layer by the separatrix of 
the resonance Hamiltonian (16). The resonance approximation (16) neglects time- 
periodic terms while the frequencies of oscillation of these terms greatly exceed the 
frequency of eigenoscillation of the resonance Hamiltonian (16) around its relevant 
elliptic point i.e. the elliptic point inside the area limited by the resonance separa- 
trix. As is well known ITSl l23l l29l ISTl l52l l55l . fast-oscillating terms acting on a 
system with a separatrix give rise to the onset of an exponentially narrow chaotic 
layer in place of the separatrix. In the present context, this means that the correction 



\aeT 

AE 



.{next) 



1 



ln(l//i)' 



(74) 
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to the maximal action I stemming from fast-oscillating corrections to the resonance 
Hamiltonian, i.e. AEf"''''^'\ is exponentially small, thus being negligible in compari- 
son with the correction Zi/s,',"^'^'^ (see (74)). 

The right wing, described in Sees. 3.2 and 3.3 above, corresponds in leading- 
order approximation to the formation of the boundary of the layer by the resonance 
trajectory tangent to the GSS curve. For the part of the right wing exponentially 
close in frequency to the frequency of the maximum, the tangent trajectory is close 
to the resonance separatrix, so that the correction stemming from fast-oscillating 
terms is exponentially small, similarly to the case of the left wing. As the frequency 
further deviates from the frequency of the maximum, the tangent trajectory further 
deviates from the resonance separatrix and the correction AE, differs from the 
exponentially small correction estimated above. It may be estimated in the following 
way. 

It follows from the second-order approximation of the averaging method fS] that 
the fast-oscillating terms lead, in the second-order approximation, to the onset of ad- 
ditional terms h^Tj{I, iff) and h^Tif,{I, iff) in the dynamic equations for slow variables 
I and iff respectively, where Tj{I, iff) and r^(/, iff) are of the order of the power-law- 
like function of l/ln(l /h) in the relevant range of I. The corresponding correction 
to the width of the chaotic layer in energy may be expressed as 



AEj"'"^= dth^Tjcoil), (75) 

where fmin and f^ax are instants corresponding to the minimum and maximum de- 
viation of the tangent trajectory from the separatrix of the unperturbed system (cf. 
Figs. 1(c) and 4(c)). The interval tmax ~ fmin may be estimated as follows: 

^max ^min ^ "i I T? (76) 

\<w>\ 

where < iff > is the value of iff averaged over the tangent trajectory. It follows from 
(16) that 



<iff>\^C0f- co{E, -8)^ ,^7; ^ (77) 



CO{Es -5) (Oq 

ln(l//2) ^ ln2(l//i) 

Taking together Eqs. (75)-(77) and allowing for the fact that 7} is of the order of 
a power-law-like function of l/ln(l//i), we conclude that 

AE^"'"'^ ^h^P{ln{l/h)), (78) 

where P{x) is some power-law-like function. 

The value Zifi/""'^ is still asymptotically smaller than the absolute value of the 
correction within the resonance approximation, \AEi"''^^\, which is of the order of 
h or /j/ln(l /h) for systems of type I or type II respectively. 

Thus, we conclude that, both for the left and right wings of the peak, (i) the cor- 
rection AEf""'^''' is determined by the correction within the resonance approximation 
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AEw"'\ and (ii) in the asymptotic limit h 0, the overall next-order correction is 
negligible in comparison with the leading term: 



AE AE AE ln(l//2) 

This estimate well agrees with results in Figs. 3-6. 



0. (79) 



3.5 Discussion 



In this section, we briefly discuss the following issues: (i) the scaled asymptotic 
shape of the peaks; (ii) peaks in the range of moderate frequencies; (iii) jumps in 
the amplitude dependence of the layer width; and (iv) chaotic transport; (v) smaller 
peaks at rational frequencies; (vi) other separatrix maps; (vii) an application to the 
onset of global chaos. 

1. Let us analyse the scaled asymptotic shape of the peaks. We consider first sys- 
tems of type I. The peaks are then described in the leading-order approximation 
exclusively within separatrix map dynamics (approximated, in turn, by the NR 
dynamics). It follows from Eqs. (32), (34), (36), (39) and (40) that most of the 
peak fir given j can be written in the universal scaled form: 

4£W(tO;) = AE^l^S ("^^-^ifOf-oiil)] , (80) 
where the universal function S{a) is strongly asymmetric: 



^y^)~\sAa) for a>0, 

It is not difficult to show that 



Si{a = 0) = 1 
d5,(a = 0) 



da 



= l-e" 



Sria 



-0) 



(82) 



da 



S{a ±oo) <^ 



a 
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Thus, the function S{a) is discontinuous at the maximum. To the left of the 
maximum, it approaches the far part of the wing (which decreases in a power- 
law-Uke way) relatively slowly while, to the right of the maximum, the function 
first drops jump-wise by a factor e and then sharply approaches the far part of the 
wing (which again decreases in a power-law-like way). 

It follows from Eqs. (80), (81), (82) and (27) that the peaks are logarithmically 
narrow, i.e. the ratio of the half-width of the peak, Aco'^-'K to COmax is logarithmi- 
cally small: 

Aa'^j^ 1 

"^^ln(8(2y-l)//.)- ^^^^ 

We emphasize that the shape (81) is not restricted to the example (14): it is valid 
for any system of type I. 

For systems of type II, contributions from the NR and from the variation of en- 
ergy during the pulse of velocity, in relation to their li dependence, are formally 
of the same order but, numerically, the latter contribution is usually much smaller 
than the former. Thus, typically, the function (81) approximates well the properly 
scaled shape of the major part of the peak for systems of type II too. 

2. The quantitative theory presented in the paper relates only to the peaks of small 
order n i.e. in the range of logarithmically small frequencies. At the same time, 
the magnitude of the peaks is still significant up to frequencies of order of one. 
This occurs because, for motion close to the separatrix, the order of magnitude of 
the Fourier coefficients remains the same up to logarithmically large numbers «. 
The shape of the peaks remains the same but their magnitude typically decreases 
(though in some cases, e.g. in case of the wave-like perturbation ll23l ISTl [52l l55l 
it may even increase in some range of frequencies). The quantitative description 
of this decrease, together with analyses of more sophisticated cases, requires a 
generalization of our theory. 

3. Apart from the frequency dependence of the layer width, our theory is also rele- 
vant to amplitude dependence: it describes the jumps ll40l in the dependence of 
the width on h and the transition between the jumps and the linear dependence. 
The values of h at which the jumps occur, h^j^lup, are determined by the same 

condition that determines a),nax in the frequency dependence of the width. The 
formulae relevant to the left wings of the peaks in the frequency dependence de- 
scribe the ranges h > h^j'J^jp while the formulae relevant to the right wings describe 

the ranges li < h^fj^^. 

4. Apart from the description of the boundaries, the approach allows us to describe 
chaotic transport within the layer In particular, it allows us to describe quantita- 
tively the effect of the stickiness of the chaotic trajectory to boundaries between 
the chaotic and regular areas of the phase space i51..52J. Moreover, the presence 
of additional (resonance) saddles should give rise to an additional slowing down 
of the transport, despite a widening of the area of the phase space involved in the 
chaotic transport. 
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5. Our approach can be generalized in order to describe smaller peaks at non-integer 
rational frequencies i.e. cOf w n/mcOr"^^ where n and m are integer numbers. 

6. Apart from Hamiltonian systems of the one and a half degrees of freedom and 
corresponding Zaslavsky separatrix maps, our approach may be useful in the 
treatment of other chaotic systems and separatrix maps (see ||29l for the most 
recent major review on various types of separatrix maps and related continuous 
chaotic systems). 

7. Finally we note that, apart from systems with a separatrix, our work may be 
relevant to nonlinear resonances in any system. If the system is perturbed by a 
weak time-periodic perturbation, then nonlinear resonances arise and their dy- 
namics is described by the model of the auxiliary time-periodically perturbed 
pendulum HO] |23l |55l |5l] |52l [U [l8l| . If the original perturbation has a single 
harmonic, then the effective perturbation of the auxiliary pendulum is necessar- 
ily a high-frequency one, and chaotic layers associated with the resonances are 
exponentially narrow ifTol l23l l55l ISTl [52l [l] [l8l while our results are irrelevant. 
But, if either the amphtude or the angle of the original perturbation is slowly 
modulated, or if there is an additional harmonic of a slightly shifted frequency, 
then the effective perturbation of the auxiliary pendulum is a low-frequency one 
B3l and the layers become much widej^ while our theoretical approach becomes 
relevant. It may allow to find optimal parameters of the perturbation for the facil- 
itation of the onset of global chaos associated with the overlap in energy between 
different-order nonlinear resonances tlOJ : the overlap may be expected to occur 
at a much smaller amplitude of perturbation in comparison with that one required 
for the overlap in case of a single-harmonic perturbation. 



4 Double-separatrix chaos 

There are many problems in physics where an unperturbed Hamiltonian model pos- 
sesses two or more separatrices. A weak perturbation of the system typically de- 
stroys the separatrices, replacing them by thin chaotic layers. As the magnitude of 
the perturbation grows, the layers become wider and, at some critical value, they 
merge with each other: this may be described as the onset of global chaos between 
the separatrices. Such a connection of regions of different separatrices is important 
for transport in the system. 

In the present section, following the paper 1431 . we consider the characteristic 
problem of the onset of global chaos between two close separatrices of a ID Hamil- 
tonian system perturbed by a time-periodic perturbation. As a characteristic example 
of a Hamiltonian system with two or more separatrices, we use a spatially periodic 
potential system with two different-height barriers per period (Fig. 7(a)): 



This should not be confused with the widening occuring with the separatrix chaotic layer in the 
original pendulum if an originally single-harmonic perturbation of a high frequency is completed 
by one more harmonic of a slightly shifted frequency: see |47| and references therein. 
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P , . , . cP — sin o) 
Ho{p,q) = Y + U{q), U{q) = ^ ^ = const <1. (84) 

This model may relate e.g. to a pendulum spinning about its vertical axis [J or 
to a classical 2D electron gas in a magnetic field spatially periodic in one of the 
in-plane dimensions P9l l50l . Interest in the latter system arose in the 1990s due 
to technological advances allowing to manufacture magnetic superlattices of high- 
quality |l6]|48l, and thus leading to a variety of interesting behaviours of the charge 
carriers in semiconductors B9l l50l l6l l48l |331 [33l . 

Figs. 7(b) and 7(c) show respectively the separatrices of the Hamiltonian sys- 
tem (1) in the p — q plane and the dependence of the frequency (O of its oscillation, 
often called its eigenfrequency, on its energy E = H(){p,q). The separatrices corre- 
spond to energies equal to the value of the potential barrier tops E^'' = (1 ^ ^')^/2 
and E^^' = (1 + 0)^/2 (Fig. 7(a)). The function co{E) possesses a local maximum 
CO,n = Co{E,„). Moreover, co{E) is close to co,„ for most of the range [E^'^E^^'] while 

sharply decreasing to zero as E approaches either Zs^^'^ or Ej^\ 

We now consider the addition of a time-periodic perturbation: as an example, 
we use an AC drive, which corresponds to a dipole 1551 |2T1 perturbation of the 
Hamiltonian: 

q = dH/dp, p = -dH/dq, (85) 
H{p,q) ^ Ho{p,q) - hqcos{(Oft). 

The conventional scenario for the onset of global chaos between the separatrices 
of the system (84)-(85) is illustrated by Fig. 8. The figure presents the evolution of 



0.8 




°0 0.2 0.4 0.6 0.8 



Energy £ 

Fig. 7 The potential [/((?), the separatrices in the phase space, and the eigenfrequency (o{E) for 
the unperturbed system (84) with <P = 0.2, in (a), (b) and (c) respectively. 
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Fig. 8 The evolution of the 
stroboscopic (at t = nlK / cOf 
with « = 0, 1 , 2, . . .) Poincare 
section of the system (84)- 
(85) with 4> = 0.2 as h grows 
while cOf = 0.3. The number 
of points in each trajectory 
is 2000. In (a) and (b), three 
characteristic trajectories are 
shown: the inner trajectory 
starts from the state } = 
{p = Q,q = ^/2} and is 
chaotic but bounded in space; 
the outer trajectory starts from 
{0(0'} = {p = 0,</ = -;r/2} 
and is chaotic and unbounded 
in coordinate; the third tra- 
jectory is an example of a 
regular trajectory separating 
the two chaotic ones. In (c), 
the chaotic trajectories mix. 

the stroboscopic Poincare section as h grows while cOf is fixed at an arbitrarily cho- 
sen value away from co„, and its harmonics. At small li, there are two thin chaotic lay- 
ers around the inner and outer separatrices of the unperturbed system. Unbounded 
chaotic transport takes place only in the outer chaotic layer i.e. in a narrow energy 
range. As li grows, so also do the layers. At some critical value ligc = hgc{(Of), the 
layers merge. This may be considered as the onset of global chaos: the whole range 
of energies between the barrier levels is involved, with unbounded chaotic transport. 
The states } = {/:> = 0, ^ = 7r/2 + 27r/} and {0^''> } = {p = 0,q = -7t/2 + 27rl} 
(where / is any integer) in the Poincare section are associated respectively with the 
inner and outer saddles of the unperturbed system, and necessarily belong to the 
inner and outer chaotic layers, respectively. Thus, the necessary and sufficient con- 
dition for global chaos onset may be formulated as the possibility for the system 
placed initially in the state {/'*^'} to pass beyond the neighbourhood of the "outer" 
states, {(9*"'} or {(9^''}, i.e. for the coordinate q to become < — 7r/2 or > 37r/2 at 
sufficiently large times t ^ 2n/(0f. 

A diagram in the h — (Of plane, based on the above criterion, is shown in Fig. 
9. The lower boundary of the shaded area represents the function hgc{(£>f). It has 
deep spikes i.e. cusp-like local minima. The most pronounced spikes are situated at 
frequencies (Of ~ (Os^'' that are slightly less than the odd multiples of (»,„, 

«i^'«a),„(2;-l), j = l,2,... (86) 

The deepest minimum occurs at coi w ©„,: the value of hgc at the minimum, = 

hgc{(Os^''), is approximately 40 times smaller than the value in the neighbouring 
pronounced local maximum of hgc{(Of) at O/^ « 1. As n increases, the «th minimum 
becomes shallower The function hgc{(Of) is very sensitive to cOf in the vicinity of 




Coordinate q h=0.1 
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Fig. 9 Diagram indicating 
the range of perturbation 
parameters (shaded) for which 
global chaos exists. The 
integration time for each point 
of the grid is 12000;^:. 




the minima: for example, a reduction of (Of from O)] w 0.4 of only 1% causes an 
increase in hgc of w 30%. 

The origin of the spikes is related to the involvement of the resonance dynamics 
in separatrix chaos, similar to that considered in Sec. 3. In particular, the minima of 
the spikes correspond to the situation when the resonances almost touch, or slightly 
overlap with, the separatrices of the unperturbed system while overlapping each 
other This is illustrated by the evolution of the Poincare section as h grows while 
(Of « coi'' (Fig. 10) and by its comparison with the corresponding evolution of 
resonance separatrices calculated in the resonance approximation (Fig. 1 1). 

Sec. 4. 1 below presents the self-consistent asymptotic theory of the minima of the 
spikes, based on an accurate analysis of the overlap of resonances with each other 
and on the matching between the separatrix map and the resonance Hamiltoinian 
(details of the matching are developed in Appendix). Sec. 4.2 presents the theory of 
the wings of the spikes Generahzations and applications are discussed in Sec. 4.3. 



4.1 Asymptotic Theory For The Minima Of The Spikes 

The eigenfrequency Co{E) stays close to its local maximum co,„ for most of the rele- 
vant range [Ej^^\Ej^^] (Fig. 7(c)). As shown below, Co{E) approaches a rectangular 
form in the asymptotic limit ^> ^ 0. Hence, if the perturbation frequency (Of is 
close to (0,„ or its odd multiples, \(0f ~ {2j — l)(Om \ ^ (Om, then the energy widths of 
nonlinear resonances become comparable to the width of the whole range between 
the barriers (i.e. E^^' — ii^^' « 2<P) at a rather small perturbation magnitude li ^ <P. 
Note that <P determines the characteristic magnitude of the perturbation required 
for the conventional overlap of the separatrix chaotic layers, when (0/ is not close 

to any odd multiple of (0,„ (Fig. 8 (c)). Thus, if o)/ ~ (Os'^\ the nonlinear resonances 
should play a crucial role in the onset of global chaos (cf. Fig. 10). 

We note that it is not entirely obvious a priori whether it is indeed possible to cal- 
culate /i'^' = hgc{(Os''^) within the resonance approximation: in fact, it is essential for 
the separatrices of the nonlinear resonances to nearly touch the barrier levels, but the 
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Fig. 10 (Color version may 
be found in the online ver- 
sion of 1431 as Fig. 5). The 
evolution of the stroboscopic 
Poincare section of the system 
(84)-(85) with <P = 0.2, as 
the amplitude h of the per- 
turbation grows, while the 
frequency remains fixed at 
(Of = 0.401. The number of 
points in each trajectory is 
2000. The chaotic trajecto- 
ries starting from the states 
{/(")} and {0<")} are drawn 
in green and blue respectively. 
The stable stationary points 
of Eq. (98) for n = 1 (i.e. 
for the Ist-order nonlinear 
resonances) are indicated by 
the red and cyan crosses. The 
chaotic layers associated with 
the resonances are indicated 
in red and cyan respectively, 
unless they merge with those 
associated with the green/blue 
chaotic trajectories. Examples 
of regular trajectories em- 
bracing the state } while 
separating various chaotic tra- 
jectories are shown in brown. 
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Fig. 11 (Color version may 
be found in the online version 
of 1431 as Fig. 6). The evolu- 
tion of the separatrices of the 
Ist-order resonances within 
the resonance approximation 
(described by (16) with n= \) 
in the plane of action / and 
slow angle ij/, for the same pa- 
rameters as in Fig. 10 (boxes 
(a), (b), (c), (d) correspond to 
those in Fig. 10). Horizontal 
levels mark the values of / 
corresponding to the barriers. 
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resonance approximation is invalid in the close vicinity of the barriers; furthermore, 
numerical calculations of resonances show that, if (Of « ft)]^', the perturbation am- 
plitude h at which the resonance separatrix touches a given energy level in the close 
vicinity of the barriers is very sensitive to (Of, apparently making the calculation of 

within the resonance approximation even more difficult. 

Nevertheless, we show below in a self-consistent manner that, in the asymptotic 
limit ^> ^ 0, the relevant boundaries of the chaotic layers lie in the range of energies 
E where (o{E) k, (o,„. Therefore, the resonant approximation is valid and it allows 
us to obtain explicit asymptotic expressions both for o]''' and hi''\ and for the wings 
of the spikes in the vicinities of o]''' . 

The asymptotic limit <P ^0 is the most interesting one from a theoretical point 
of view because it leads to the strongest facilitation of the onset of global chaos, 
and it is most accurately described by the self-contained theory. Most of the theory 
presented below assumes this limit and concentrates therefore on the results to the 
lowest (i.e. leading) order in the small parameter. 

On the applications side, the range of moderately small <P is more interesting, 
since the chaos facilitation is still pronounced (and still described by the asymptotic 
theory) while the area of chaos between the separatrices is not too small (compa- 
rable with the area inside the inner separatrix): cf. Figs. 7, 8 and 10. To increase 
the accuracy of the theoretical description in this range, we estimate the next-order 
corrections and develop an efficient numerical procedure allowing for further cor- 
rections. 



4,1.1 Resonant Hamiltonian and related quantities 

Let (Of be close to the «th odc0 harmonic of a)„,, « = (2y — 1). Over most of 

the range [Ej^^\Ej^^], except in the close vicinities of E^^' and Ej^\ the «th har- 
monic of the eigenoscillation is nearly resonant with the perturbation. Due to 
this, the (slow) dynamics of the action / = I{E) = {2k)^^ § dqp and the angle ^ 
ifTOl l23l l55l ISTl l52l |4T1 l2Tll can be described by means of a resonance Hamiltonian 
similar in form to (16). The lower integration limit in the expression for H may be 
chosen arbitrarily, and it will be convenient for us to use presently /(£■,„) (instead of 
I{Es) in (16)) where E„, is the energy of the local maximum of (o{E) (Fig. 7(c)). To 
avoid confusion, we write the resonance Hamiltonian explicitly below after making 
this change; 

H{I,\fr)= [ dl{n(0-(0f) - nliqnCos{\ff) (87) 

J HE,,, ) 

= n{E - E,„) - (Of(I - I{E„,)) - nhq„cos{\ff) , 



Even harmonics are absent in the eigenoscillation due to the symmetry of the potential. 
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rE f\£ 

I = I{E)= E = HQ{p,q), 

iff = nyz-cOft, 

fi da 
\(/=7Z + sign{p)(o{E) =+2nl, 

A™„(£) ^2{E-U[q)) 

2 f^l^ 

q„ = q„(E)^- d\i/q{E,\l/)cos{n\l/), 



K Jo 

|«a)- O/l <c 0), « = 2y-l, y = 1,2,3,... 

Let us derive explicit expressions for various quantities in (87). In the unper- 
turbed case (h = 0), the equations of motion (85) with Hq (84) can be integrated 
ll50ll (see also Eq. (144) below), so that we can find (o{E): 



where 



, , 7r(2£)'/4 , 1 (V2l+l)2-02 

'^(^) = ^Z[r' '=2V — — ' ^''^ 

d0 



is the complete elliptic integral of first order |2|. Using its asymptotic expression, 



2 \l-k^ 
we derive Co{E) in the asymptotic limit ^> 0: 



co{E) ~ —. ^- ^, (89) 



In 



64 



[0-AE){>l>+AE) 



AE = E-^, \AE\<0, 0^0. 

As mentioned above,the function Co{E) (89) remains close to its maximum 

(0,n= max -J— (90) 

[4",£f'] 21n(8/0) 

(I 2) 

for most of the interbarrier range of energies [1/2— c?>,l/2 + ^>] (note that is^ w 
1/2 =F ^> to first order in ^>.); on the other hand, in the close vicinity of the barri- 
ers, where either |ln(l/(l -AE/<P))\ or |ln(l/(l +AE/<P))\ become comparable 
with, or larger than, ln(8/^'), Co{E) decreases rapidly to zero as \AE\ 0. The 
range where this takes place is ~ 0^, and its ratio to the whole interbarrier range, 
20, is ~ i.e. it goes to zero in the asymptotic limit ^ 0:in other words, (o{E) 
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approaches a rectangular form. As it will be clear from the following, it is this al- 
most rectangular form of (o{E) which determines many of the characteristic features 
of the global chaos onset in systems with two or more separatrices. 

One more quantity which strongly affects {(OsJis) is the Fourier harmonic q„ = 
q,i{E). The system stays most of the time very close to one of the barriers. Consider 
the motion within one of the periods of the potential U{q), between neighboring 

upper barriers [q[2,q[fij] where ^J^' = q[2 + 2n. If the energy E = 1/2 + AE lies 

in the relevant range [Ejy^\Ej^'^], then the system will stay close to the lower barrier 

qih = ql^^ + n for a tim^B 

during each period of eigenoscillation, while it will stay close to one of the upper 

fi 

barriers ^j,^" = qii,±7t for most of the remainder of the eigenoscillation. 

Hence, the function q{E,\i/) — qn, may be approximated by the following piecewise 
even periodic function: 



n at ^e[o,^j2M^]u[n-^j^. 
at ^e]i^,n-§jli 



qiE,w)^q,b^{^ ::i.'rvj ' (93) 



q{E,-\i/)-qih = q{E,w)-qih, q{E,w±^^i) = q{E,w), / = 1,2,3,... 

Substituting the above approximation for q{E, y/) into the definition of q^ (87), one 
can obtain: 



qij-i = q2j-i{E) = - ^ . sm 



(2;-l);r/2 



(94) 



27-1 

0^0, ^2,/ = 0, j= 1,2,3,... 
At barrier energies, q^j- \ takes the values 

92;-i(4'^)=0, q^._,[Ef^) = -[-iy^-l—^. (95) 

As E varies in between its values at the barriers, ^2j-i varies monotonically if 
7 = 1 and non-monotonically otherwise (cf. Fig. 16). But in any case, the significant 



We omit corrections ~ (ln(l/^>)) ' here and in Eq. (92) since they vanish in the asymptotic 
limit "J) ^ 0. 
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variations occur mostly in the close vicinity of the barrier energies E^'^ and E^^'' 

while, for most of the range [e\P ^E^'"'], the argument of the sine in Eq. (94) is close 
to n/A and ^27-1 is then almost constant: 



\+AE/^ 



2ti 



In 



\-AE/^ 



V2 
<21n 



= 1,2,3, 



(96) 



where [. . .] means the integer part. 

In the asymptotic limit <P -^0, the range of AE for which the approximate equal- 
ity (96) for qij-i is valid approaches the whole range ] — 

We emphasize that \qn\ determines the "strength" of the nonUnear resonances: 
therefore, apart from the nearly rectangular form of Co{E), the non-smallness of |^„| 
is an important additional factor strongly facilitating the onset of global chaos. 

We shall need also an asymptotic expression for the action /. Substituting Co{E) 
(89) into the definition of I{E) (87) and carrying out the integration, we obtain 



AE\n 



1(1 12)- 



646-^ 



(97) 



4.1.2 Reconnection of resonance separatrices 



We now turn to analysis of the phase space of the resonance Hamiltonian (87). The 
evolution of the Poincare section (Fig. 10) suggests that we need to find a separatrix 
of (87) that undergoes the following evolution as h grows: for sufficiently small h, 
the separatrix does not overlap chaotic layers associated with the barriers while, for 
h > li^c{(Of), it does overlap them. The relevance of such a condition will be further 
justified. 

Consider cOf ~ n(0,„ with a given odd n. For the sake of convenience, let us write 
down the equations of motion (87) explicitly: 



i = = -nhqnsm{\ir), \j/ = — — = - ©/■ - -^cos(v>'). (98) 
a^i ol dl 

Any separatrix necessarily includes one or more unstable stationary points. The sys- 
tem of dynamic equations (98) may have several stationary points per In interval of 
iff. Let us first exclude those points which are irrelevant to a separatrix undergoing 
the evolution described above. 

Given that ^„(£'^'') = 0, there are two unstable stationary points with / corre- 
sponding to £■ = /i^'' and xjf = ±7z/2. They are irrelevant because, even for an 
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Fig. 12 (Color version may 
be found in the online ver- 
sion of 1431 as Fig. 7). A 
schematic example illustrat- 
ing the graphical solutions 
of Eqs. (99) for o = 1, as in- 
tersections of the curve ft)(/) 
(thick solid red line) with 
the curves CO/ ±Mt/„(/)/d/ 
(thin solid green lines). The 
solutions con'esponding to 
the lower and upper relevant 
saddles (defined by Eq. (100)) 
are marked by dots and by 
the labels si and su respec- 
tively (we do not mark other 
solutions because they are 
irrelevant). 

infinitely small h, each of them necessarily lies inside the corresponding barrier 
chaotic layer 

If £■ 7^ E^i^\ then ^„ 7^ 0, so / = only if \J/ is equal either to or to n. Substituting 
these values into the second equation of (98) and putting \J/ = 0, we obtain the 
equations for the corresponding actions: 

X^{I)=n(0-(OfTnhdq„/dJ = Q, (99) 

where the signs "-" and "+" correspond to ij/ = and iff ~ n respectively. A typical 
example of the graphical solution of equations (99) for « = 1 is shown in Fig. 12. 
Two of the roots corresponding to xfr ~ n are very close to the barrier values of / 
(recall that the relevant values of h are small). These roots arise due to the divergence 
of dq/dl as / approaches any of the barrier values. The lower/upper root corresponds 
to a stable/unstable point, respectively. However, for any «, both these points and the 
separatrix generated by the unstable point necessarily lie in the ranges covered by 
the barrier chaotic layers. Therefore, they are also irrelevant For n > 1 , the number 
of roots of (99) in the vicinity of the barriers may be larger (due to oscillations of the 
modulus and sign of dq„/dl in the vicinity of the barriers) but they all are irrelevant 
for the same reason, at least to leading-order terms in the expressions for the spikes' 
minima. 

Consider the stationary points corresponding to the remaining four roots of equa- 
tions (99). Just these points are conventionally associated with nonlinear resonances 
|[l0l|23l|55l|5l]|52l|4t|. It follows from the analysis of equations (98) linearized near 
the stationary points (cf. |[l0l|23l|55l|5ll|52llSl), two of them are stable (elliptic) 




' ' For sufficiently small (p and li, the separatrix generated by the unstable point forms the boundary 
of the upper chaotic layer, but this affects only the higher-order terms in the expressions for the 
spikes miniina (see below). 
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point0, while two others are unstable (hyperbolic) points, often called saddles. 
These saddles are of central interest in the context of our work. They belong to 
the separatrices dividing the I — iff plane for regions with topologically different 
trajectories. 

We shall identify the relevant saddles as those with the lower action/energy (using 
the subscript "iZ") and upper action/energy (using the subscript "sm"). The positions 
of the saddles in the /— ij/ plane are defined by the following equations (cf. Figs. 1 1 
and 12); 



g = sgn(^„ (/,„,,,)) = sgn ((-!)[?]) , (100) 

CU.v/ d/.s„ 

where [...] means an integer part, X±{I) are defined in Eq. (99) while 7^/ and 7™ are 
closer to I{E,n) than any other solution of (100) (if any) from below and from above, 
respectively. 

Given that the values of h relevant to the minima of the spikes asymptotically 
approach in the asymptotic limit 4" ^ 0, one may neglect the last term in the def- 
inition of Xzf^ in Eq. (99) in the lowest-order approximatiorQ so that the equations 

— reduce to the simplified resonance condition 

n(0{lsu.si) = «/• (101) 

Substituting here Eq. (89) for ©, we obtain explicit expressions for the energies in 
the saddles: 



i±4£(i), (102) 

AE'^^'^ = y^$>2-64exp^-^^, (Of < n(Om. 

The corresponding actions 7„, are expressed via Esu,si by means of Eq. (97). 

For (Of ~ n(0,n, the values of E^u,si (102) lie in the range where the expression 
(96) for q„ holds true. This will be confirmed by the results of calculations based on 
this assumption. 

Using (100) for the angles and (102) for the energies, and the asymptotic expres- 
sions (89), (96) and (97) for (o{E), qn{E) and 7(7:) respectively, and allowing for 



In the Poincare sections shown in Fig. 10, the points which correspond to such stable points of 
equations (98) are indicated by the crosses. 

As will become clear in what follows, the remaining terms are much larger in the asymptotic 
limit than the neglected term: cf. the standard theory of the nonlinear resonance 11011231151115211551 . 



A New Approach To The Treatment Of Separatrix Chaos And Its Applications 



43 



the resonance condition (101), we obtain explicit expressions for the values of the 
Hamiltonian (87) at the saddles: 

+ hV2. (103) 

As the analysis of simulations suggests and as it is self-consistently shown fur- 
ther, one of the main conditions which should be satisfied in the spikes is the overlap 
in phase space between the separatrices of the nonlinear resonances, which is known 
as separatrix reconnection PTl [191 |20l [TTl [T2l |26l . Given that the Hamiltonian H is 
constant along any trajectory of the system (87), the values of H in the lower and 
upper saddles of the reconnected separatrices are equal to each other: 

H,,=H,„. (104) 

This may be considered as the necessary and sufficient condition for the recon- 
nection. Taking into account that Hsi = —Hsu (see (103)), it follows from (104) that 

^,/=i?,„=0. (105) 
Explicitly, the relations in (105) reduce to 



-nsii — 

7t 



2zi£(i)-01n 



+ 



h = h{(Of) 

n = 1,3,5, 



_COf_ 

V2% 



01n I 



+ 



Jl TT 

-64exp( ), 



< ft)„, -(Of/n<ii, 0}„, 



(106) 



21n(8/0)' 



The function li{(Of) (106) decreases monotonically to zero as cOf grows from to 
n(Om, where the line abruptly stops. Fig. 15 shows the portions of the lines (106) 
relevant to the left wings of the 1st and 2nd spikes (for = 0.2). 



4.1.3 Barrier chaotic layers 

The next step is to find the minimum value of h for which the resonance separatrix 
overlaps the chaotic layer related to a potential barrier. With this aim, we study 
how the relevant outer boundary of the chaotic layer behaves as li and cOf vary. 
Assume that the relevant (Of is close to ncOm while the relevant li is sufficiently 
large for (o{E) to be close to o,,, at all points of the outer boundary of the layer (the 
results will confirm these assumptions). Then the motion along the regular trajectory 



Eq. (104) is the sufficient (rather than just necessary) condition for separatrix reconnection 
since there is no any other separatrix which would lie in between the separatrices generated by the 
saddles "si" and "su". 
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infinitesimally close to the layer boundary may be described within the resonance 
approximation (87). Hence the boundary may also be described as a trajectory of 
the resonant Hamiltonian (87). This is explicitly proved in the Appendix, using a 
separatrix map analysis allowing for the validity of the relation Co{E) w ©„, for all E 
relevant to the boundary of the chaotic layer The main results are presented below. 
For the sake of clarity, we present them for each layer separately, although they are 
similar in practice. 

4. 1 .3. 1 Lower Layer 

Let (Of be close to any of the spikes' minima. 

One of the key roles in the formation of the upper boundary of the layer is played 
by the angle-dependent quantity 5/ 1 sin( | which we call the generalized separatrix 
split (GSS) for the lower layer, alluding to the conventional separatrix split ISTI for 
the lower layer 5/ = |e('""') (a)/)|/! with e'''™' given by Eq. (17213 (cf. also (4)). 
Accordingly, we use the term "lower GSS curve" for the following curve in the 
I — \fr plane: 



4.1.3.1.1 Relatively Small Values Ofh 

If h < he) {(Of), where the critical value hi!) {(Of) is determined by Eq. (125) (its 
origin will be explained further), then there are differences in the boundary for- 
mation for the frequency ranges of odd and even spikes. We describe these ranges 
separately. 



In this case, the boundary is formed by the trajectory of the Hamiltonian (87) 
tangent to the GSS curve (see Fig. 22(a); cf. also Figs. 1(c), 13(a), 14(b), 14(c)). 
There are two tangencies in the angular range ] — ;r,;r[: they occur at the angles 
±1///'^ where yfj'^ is determined by Eq. (182). 

In the ranges of h and (0/ relevant to the spike minimum, the asymptotic expres- 
sions for 5; and are; 



/ = /«s(V/)^/(4') + 5,|sin(vA)|). 



(107) 



1. Odd spikes 



5i = Vlnh, 



(108) 



(109) 



81n(l/0) 



2 



The quantity 5; may also be interpreted as the magnitude of the con'esponding Melnikov integral 
I10ll23ll55ll51ll52l . sometimes called as the Poincare-Melnikov integral 1291 . 
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Hence, the asymptotic value for the deviation of the tangency energy from the 
lower barrier reduces to: 

E^'^-Ei^^ ^5,sinifP)^^-^ '] . (110) 

2 y/ln{l/<P)/n 

The minimum energy on the boundary, corresponds to xfr ~Qoi n for even 
or odd values of [n/4] respectively. Thus, it can be found from the equality 

^(^('El!]n),V^-^(l-(-l)['b/2) (Ill) 

At 0, Eq. (Ill) yields the following expression for the minimal deviation 
of energy on the boundary from the barrier: 

.(/) ^ (/) _ (1) _ . (/) (112) 

In the context of the onset of global chaos, the most important property of the 
boundary is that the maximum deviation of its energy from the barrier, 5mix, should 
greatly exceed both d'fl^ and 5/. As h li^r , the maximum of the boundary ap- 
proaches the saddle "si". 

2. Even spikes 

In this case, the Hamiltonian (87) possesses saddles 'V in the close vicinity to 
the lower barrier (see Fig. 22(b)). Their angles differ by n from those of "si": 

V/v = 7r ^-^-^ V2nm, m = 0,±1,±2, . . . , (113) 

while the deviation of their energies from the barrier still lies in the relevant (reso- 
nant) range and reads, in the lowest-order approximation, 

71 h 

^■^■"2x/2ln(ln(l/0))- ^^^^^ 

The lower whiskers of the separatrix generated by these saddles intersect the 
GSS curve while the upper whiskers in the asymptotic limit do not intersect it (Fig. 
22(b)). Thus, it is the upper whiskers of the separatrix which form the boundary 
of the chaotic layer in the asymptotic limit. The energy on the boundary takes the 
minimal value right on the saddle "s", so that 

^ (^^^^ 

™" 2x/2ln(ln(l/0))" ^ ' 
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(b) 



Fig. 13 (Color version may be found in tlie online version of (43] as Fig. 8). (a) Chaotic layers 
(shaded in green and blue, for the upper and lower layers respectively) in the plane of action / and 
slow angle ij/, as described by our theory. Parameters are the same as in Figs. 10(b) and 11(b). 
The lower and upper boundaries of the figure box coincide with and respectively. 

The resonance separatrices are drawn by the cyan and red solid lines (for the lower and upper 
resonances respectively). Dashed green and blue lines mark the curves / = /^jg (t//) = /(£ = E^p + 



5/1 sin(ij/)|) and / = = ~ ^b^' ~ ^'^(v/)!) respectively, where 5/ and 5„ are the 

values of the separatrix split related to the lower and upper barrier respectively. The upper boundary 
of the lower layer is formed by the trajectory of the resonant Hamiltonian system (87) tangent to 
the curve / = Iq\s{y)- The lower boundary of the upper layer is formed by the lower part of the 
upper (red) resonance separatrix. The periodic closed loops (solid blue lines) are the trajectories 
of the system (87) tangent to the curve /gss(V'): they form the boundaries of the major stability 
islands inside the upper chaotic layer, (b) Comparison of the chaotic layers obtained from computer 
simulations (dots) with the theoretically calculated boundaries (solid lines) shown in the box (a). 



7(2) . 



Similar to the case of the odd spikes, the maximal deviation of the energy from 
the barrier (measured along the boundary) greatly exceeds both dj^l^ and 5/. As 
h^h['), the maximum of the boundary approaches the saddle "sf. 

4.1.3.1.2 Relatively Large Values Ofh 

\i h> he) {(Of), the previously described trajectory (either the tangent one or the 
separatrix, for the odd or even spike ranges respectively) is encompassed by the 
separatrix of the lower nonlinear resonance and typically forms the boundary of a 
major stability island inside the lower layer (reproduced periodically in iff with the 
period In). The upper outer boundary of the layer is formed by the upper part of 
the resonance separatrix. This may be interpreted as the absorption of the lower 
resonance by the lower chaotic layer. 



4.1.3.2 Upper Layer 

Let (Of be close to any of the spikes' minima. 

One of the key roles in the formation of the lower boundary of the layer is played 
by the angle-dependent quantity 5^1 sin(vA)| which we call the generalized separa- 
trix split (GSS) for the upper layer; 5„ is the separatrix split for the upper layer: 
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5u — \e^"''\(Of)\h with £("p) given by Eq. (204). Accordingly, we use the term "up- 
per GSS curve" for the following curve in the / — i// plane: 

I = 4s(r) ^ ' - sin(VA)l). (116) 



4.1.3.2.1 Relatively Small Values Ofh 

If h < hi"\(Of), where the critical value hcr\(Of) is determined by Eq. (126) (its 
origin will be explained further), then there are some differences in the boundary 
formation in the frequency ranges of odd and even spikes: for odd spikes, the for- 
mation is similar to the one for even spikes in the lower-layer case and vice versa. 

1. Odd spikes 

In the case of odd spikes, the Hamiltonian (87) possesses saddles "i" in the close 
vicinity of the upper barrier, analogous to the saddles "i" near the lower barrier in 
the case of even spikes. Their angles are shifted by n from those of "i": 

xj/s = + 7Z = 71 \-27Zm, m = 0,±l,±2,... (117) 

The deviation of their energies from the upper barrier coincides, in the lowest-order 
approximation, with 5s: 

^■^=^^-^^ln(ln(l/0))- ^'''^ 

The upper whiskers of the separatrix generated by these saddles intersect the 
upper GSS curve while the lower whiskers in the asymptotic limit do not intersect 
it. Thus, it is the lower whiskers of the separatrix which form the boundary of the 
chaotic layer in the asymptotic limit. The deviation of the energy from the upper 
barrier takes its minimal value (measured along the boundary) right on the saddle 



S^">=S~=— - (119) 

™" 2V2ln(ln(l/0))' ^ ' 

The maximal deviation of the energy from the barrier (along the boundary) 
greatly exceeds both 5^^^ and 5„. As h h["), the maximum of the boundary ap- 
proaches the saddle "su". 

2. Even spikes 

The boundary is formed by the trajectory of the Hamiltonian (87) tangent to the 
GSS curve. There are two tangencies in the angle range ] — n, 7z[: they occur at the 
angles ii///"' where i//"-' is determined by Eq. (202). 

In the ranges of h and (Of relevant to the spike minimum, the expressions for 5„ 

and XfTf""^ in the asymptotic Umit are similar to the analogous quantities in the 
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lower-layer case: 



5u = \/2nh, 



~ ill 
W 



^ : 



81"^) 



(120) 
(121) 



Hence, the asymptotic value for the deviation of the tangency energy £, from the 
upper barrier reduces to: 



e;(2)_ p(«)o 
h t "u 



7r3/2 h 



2 Vln(l/^)/« 



(122) 



The maximal energy on the boundary, E^lx, corresponds to \J/ = tt ( 1 + ( — 1 ) ) /2. 
Thus, it can be found from the equality 



H(/ = /(£;n7x),V/-;r(l + (-l)'''/'V2) (123) 



At <J> ^ 0, Eq. (123) yields the following expression for the minimal deviation 
of energy from the barrier (measured along the boundary): 

(2)_ H (2) ^^^^ (124) 



4.1.3.2.2 Relatively Large Values Ofh 

If h > her {(Of) (cf. Fig. 13(a)), the previously described trajectory (either the tan- 
gent one or the separatrix, for the even and odd spikes ranges respectively) is en- 
compassed by the separatrix of the upper nonlinear resonance and typically forms 
the boundary of a major stability island inside the upper layer (reproduced periodi- 
cally in ij/ with the period 2k). The lower outer boundary of the layer is formed in 
this case by the lower part of the resonance separatrix. This may be interpreted as 
the absorption of the upper resonance by the upper chaotic layer 



The self-consistent description of chaotic layers given above, and in more detail 
in the Appendix, is the first main result of this section. It provides a rigorous basis 
for our intuitive assumption that the minimal value of h at which the layers overlap 
corresponds to the reconnection of the nonlinear resonances with each other while 
the reconnected resonances touch one of the layers and also touch/overlap another 
layer It is gratifying that we have obtained a quantitative theoretical description of 
the chaotic layer boundaries in the phase space, including even the major stabil- 
ity islands, that agrees well with the results of simulations (see Fig. 13(b)). To the 
best of our knowledge it was the first ever B3l quantitative description of the layer 
boundaries in the phase space. 
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4.1.4 Onset of global chaos: the spikes' minima 

The condition for the merger of the lower resonance and the lower chaotic layer may 
be written as 

H{I = I{E = eI'^ + 5^J, r = ^(1 - (-1)['""V2) = H,. (125) 

The condition for the merger of the upper resonance and the upper chaotic layer 
may be written as 

H{I = I{E = £f) - 8^,Y = 7r(l + {-lt^'^)/2) = //,„. (126) 

For the onset of global chaos related to the spike minimum, either of Eqs. (125) 
and (126) should be combined with the condition of the separatrix reconnection 
(104). Let us seek first only the leading terms of hg = hs{'P) and (Os = (Os{'P). Then 
(104) may be replaced by its lowest-order approximation (105) or, equivalently, 
(106). Using also the lowest-order approximation for the barriers (Ej^ ' ' w 1 /2=p0), 
we reduce Eqs. (125), (126) respectively to 

H{I = I{E = 1/2-0 + dZ), r = ?r(l - (-l)["/41)/2) = 0, (127) 

H{I = I{E = 1/2 + 0- w=n{l + {-l)^"''^)/2) = 0, (128) 

where Sj^j^^ is given by (1 12) or (1 15) for the odd or even spikes respectively, while 
is given by (1 19) or (124) for the odd or even spikes respectively. 
To the leading order, the solution (/il'', oi'' ) of the system of equations (106),(127) 
and the solution {h["\coi"^) of the system of equations (106), (128) turn out identi- 
cal. For the sake of brevity, we derive below just {hi'\coP), denoting the latter, in 
short, as (/is, £0,0 

The system of algebraic equations (106) and (127) is still too complicated for us 
to find its exact solution. However, we need only the lowest-order solution - and this 
simplifies the problem. Still, even this simplified problem is not trivial, both because 
the function hs{0) turns out to be non-analytic and because Zifi''' in (106) is very 
sensitive to (Of in the relevant range. Despite these difficulties, we have found the 
solution in a self-consistent way, as briefly described below. 

Assume that the asymptotic dependence hs{'P) is; 



With account talcen of the next-order corrections, the spike minimum (/fj.fOj) coincides with 
(It'P , oi'' ) in case of an odd spike, or with {h["^ , co'"' ) in case of an even spike. This occurs because, 
in case of odd spikes, \g„{E)\ increases/decreases as E approaches the relevant vicinity of the 
upper/lower barrier, while it is vice versa in the case of even spikes. And the larger |f/„| the further 
the resonance separatrix extends: in other words, the reconnection of the barrier chaotic layer with 
the resonance separatrix requires a smaller value of h at the baiTier where \qn\, in the relevant 
vicinity of the barrier, is larger 
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hs = a ^ (129) 
ln(4e/^>) 

where the constant a may be found from the asymptotic solution of (106), (127), 
(129). 

Substituting the energies £ = 1 /2 - ^ + and £ = 1 /2 + - S^ax in (89) 
and taking account of (1 12), (1 15), (1 19), (124) and (129), we find that, both for the 
odd and even spikes, the inequality 

»„, -«(£') <C fi)„, (130) 
holds over the whole relevant range of energies, i.e. for 

+ (131) 

Thus, the resonant approximation is valid over the whole range (131). Eq. (96) for 
qn{E) is valid over the whole relevant range of energies too. 

Consider Eq. (127) in an explicit form. Namely, we express (Of from (127), using 
Eqs. (87), (96), and (97), using also (1 12) or (115) for odd or even spikes and (129): 

^ 21n(|)\ W{^^/<P))] 

We emphasize that the value of enters explicitly only the term 0{. . .) while, as 
is clear from the consideration below, this term does not affect the leading terms in 
(lis, (Os). Thus, SjIj'jj, does not affect the leading term of CO., at all, while it affects the 
leading term of only implicitly: Sj^^^ lies in the range of energies where nq„{E) « 
v^. This latter quantity is present in the second term in the curly brackets in (132) 
and, as becomes clear from further consideration, lis is proportional to it. 

Substituting (132) into the expression for AE'- ^^ in (106), using (129) and keeping 
only the leading terms, we obtain 

2V2 



4£'') = 0v/l-4e'-2, c=— fl. (133) 

n 

Substituting AE'^^^ from (133) into Eq. (106) for h{cOf) and allowing for (129) once 
again, we arrive at a transcendental equation for c: 

lnfT^^7^)-2;KW=c. Vl-4e'-2. (134) 



,1-Z(c) 

An approximate numerical solution of Eq. (134) is: 

c~ 0.179. (135) 

Thus, the final leading-order asymptotic formulae for (0/ and li in the minima of 
the spikes are the following: 
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«=1,3,5,..., 0^0, 

where the constant c ~ 0.179 is the solution of Eq. (134). 

The self-consistent derivation of the explicit asymptotic formulae for the minima 
of hgc{(£>f) constitutes the second main result of this section. These formulae allow 
one to predict immediately the parameters for the weakest perturbation that may 
lead to global chaos. 



4.1.5 Numerical example and next-order corrections 

For = 0.2, the numerical simulations give the following values for the frequencies 
at the minima of the first two spikes (see Fig. 9): 

coi'' w 0.4005 ±0.0005, w 1.24 ±0.005. (137) 

By the lowest-order formula (136), the values are: 

a)fo'w0.393, a)]o'wl.l8, (138) 

in rather good agreement with the simulations. 

The next-order correction for Oy can immediately be found from Eq. (132) for 
(Of and from Eq. (136) for hso, so that 



c 



«7r ( 1 ' 



21n(t) ' "=^'^'^'- ''''' 

The formula (139) agrees with the simulations even better than the lowest-order 
approximation: 

cofj' « 0.402, (O-J^ « 1 .21 . (140) 
For h in the spikes minima, the simulations give the following values (see Fig. 

9): 

w 0.0049, w 0.03. (141) 

The values according to the lowest-order formula (52) are: 

« 0.0032, « 0.01. (142) 

The theoretical value li'^ gives reasonable agreement with the simulation value ' . 

(2) 

The theoretical value gives the correct order of magnitude for the simulation 
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(2) 

value hs . Thus, the accuracy of the lowest-order formula (136) for lis is much lower 
than that for (Os'. this is due to the steepness of hgc{(Of) in the ranges of the spikes 
(the steepness, in turn, is due to the flatness of the function Co{E) nearits maximum). 
Moreover, as the number of the spike j increases, the accuracy of the lowest-order 
value h^so significantly decreases. The latter can be explained as follows. For the 
next-order correction to h[.Q , the dependence on <P reads as: 



si "sO „ ^ 



hiJ) ln(4e/0)- 



(143) 



At least some of the terms in this correction are positive and proportional to h[.Q (e.g. 
due to the difference between the exact equation (99) and its approximate version 
(101)), while h^:'^ is proportional to « = 2j — 1. Thus, for <P ~ 0.2, the relative 
correction for the 1st spike is ~ 0.25 while the correction for the 2nd spike is a 
few times larger i.e. ^ 1. But the latter means that, for <P — 0.2, the asymptotic 

(2) 

theory for the 2nd spike cannot pretend to be a quantitative description of h] , but 
only provides the correct order of magnitude. Besides, if n > 1 while <P exceeds 
some critical value, then the search of the minimum involves Eq. (150) rather than 
Eq. (104), as explained below in Sec. 4.2 (cf. Figs. 15(b) and 16). Altogether, this 
explains why lii^^ is larger than h[^^ only by 50% while /ip-* is larger than /i'q^ by 
200%. 

To provide a consistent explicit derivation of the correction to h[:Q is compli- 
cated. A reasonable alternative may be a proper numerical solution of the algebraic 
system of Eqs. (10413 and (125) for the odd spikes, or (126) for the even spikes . 
To this end, in Eqs. (104)'^ and (125) or (126) we use: (i) the exact values of the 
saddle energies obtained from the exact relations (100) instead of the approximate 
relations (101); (ii) the exact formula (88) for Co{E) instead of the asymptotic ex- 
pression (89); (iii) the exact functions (E) instead of the asymptotic formula (86); 
(iv) the relation (111) and the calculation of the "tangent" state (i///'',//'') by Eqs. 
(172), (183) for the odd spikes, or relation (123) and the calculation of the "tangent" 
state (i/a/"',//"') by Eqs. (202)-(204) for the even spikes. Note that, to find the ex- 
act function q„{E), we substitute into the definition of q„{E) in (87) the expliciQ 
solution for q{E,\i/): 



q{E,\jr) = arcsin 



' ri-V2E + <P 



1-77 

q{E,\l/)^7:~q{E,n-\l/) 




For n > 1, it is also necessary to check if the solution lies above the line (150). If it does not, 
then (104) should be replaced here by (150). 

In the general case of an arbitrary potential U{q), when the explicit expression for q{E, \(/) and 
Co{E) cannot be obtained, these functions can be calculated numerically. 
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q{E,y/) =q{E,2n-\ir) for \i/e[n,2n], 

ri = - {V2E - cp + l)sn- i —xj/j , (144) 

where sn(;c) is the elliptic sine 121 with the same modulus k as the full elliptic integral 
K defined in (88). The numerical solution described above gives: 



: 0.401, (/4''') « 0.005 



uwn \ / num 



\ / fiiim V / rtuin 



(145) 

; 0.052 . 



The agreement with the simulation results is: (i) excellent for cOs for the both 
spikes and for for the 1st spike, (ii) reasonable for /i, for the 2nd spike. Thus, 
if is moderately small, a much more accurate prediction for hs than that by the 
lowest-order formula is provided by the numerical procedure described above. 



4.2 Theory of the Spikes' Wings 

The goal of this section is to find the mechanisms responsible for the formation of 
the spikes' wings (i.e. the function Iigc{c0f) in the ranges of cOf slightly deviating 

from 0)]'''), and to provide for their theoretical description. 

Before developing the theory, we briefly analyze the simulation data (Fig. 9), 
concentrating on the 1st spike. The left wing of the spike is smooth and nearly 
straight. The initial part of the right wing is also nearly straighQ though less steep. 
But at some small distance from o]'' its slope changes jump-wise by a few times: 
compare the derivative dlii,^/ dcOf w 0. 1 at = 0.4 ~ 0.41 (see the left inset in Fig. 
9) and dhgc/dcOf ~ 0.4 at cO/ = 0.45 ~ 0.55 (see the main part of Fig. 9). Thus, even 
prior to the theoretical analysis, one may assume that there are a number of different 
mechanisms responsible for formation of the wings. 

Consider the arbitrary jth spike. We have shown in the previous section that, in 
the asymptotic hmit ^> ^ 0, the minimum of the spike corresponds to the intersec- 
tion between the line (104) with (125) or (126) for odd or even spikes respectively. 
We recall that: (i) Eq. (104) corresponds to the overlap in phase space between non- 
linear resonances of the same order n = 2y — 1; (ii) Eq. (125) or (126) corresponds 
to the onset of the overlap between the resonance separatrix associated respectively 
with the lower or upper saddle and the chaotic layer associated with the lower or 
upper potential barrier; (iii) for cOf = COs'\ the condition (125) or (126) also guaran- 
tees the overlap between the upper or lower resonance separatrix, respectively, and 
the chaotic layer associated with the upper or lower barrier'^. 



Provided hgc{(£>f) is smoothed over small fluctuations. 
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If (Of becomes slightly smaller than (Of , the resonances shift closer to the bar- 
riers while moving apart from each other Hence, as h increases, the overlap of the 
resonances with the chaotic layers associated with the barriers occurs earlier than 

with each other Therefore, at < oi^' — (Of <^ (0,n, the function hgc{(Of) should 
correspond approximately to the reconnection of resonances of the order « = 2y — 1 
(Fig. 14(a)). Fig. 15(a) demonstrates that even the asymptotic formula (106) for the 
separatrix reconnection line fits the left wing of the 1st spike quite well, and that the 
numerically calculated line (104) agrees with the simulations perfectly. 

If (Of becomes slightly larger than OJ.P^ then, on the contrary, the resonances shift 
closer to each other and further from the barriers. Therefore, the mutual overlap of 
the resonances occurs at smaller h than the overlap between any of them and the 
chaotic layer associated with the lower/upper barrier (cf. Figs. 10(c) and 10(d) as 
well as 11(c) and 11(d)). Hence, it is the latter overlap which determines the func- 
tion hgc{(Of) in the relevant range of (Of (Fig. 14(b)). Fig. 15 shows that hgc{(Of) is 

indeed well-approximated in the close vicinity to the right from oi''' by the numer- 
ical solution of Eq. (125) or (126), for an odd or even spike respectively and, for the 
1st spike and the given 4>, even by its asymptotic form. 



h = h{(Of) = n- 



{21n(t)+ln(|^)}-24£<^' 



2V2 



= .l^'--64sxp{- — ),n = 2j-l, |cO/-(oi^''| < (0,n 



(Of 



(146) 



The mechanism described above determines hgc{(Of) only in the close vicinity 

of (Os^\ If (Of/ti becomes too close to (0,„ or exceeds it, then the resonances are 
not of immediate relevance: they may even disappear or, if they still exist, their 
closed loops shrink, so that they can no longer provide for connection of the chaotic 
layers in the relevant range of h. At the same time, the closeness of the frequency 




Fig. 14 (Color version may be found in the online version of [43] as Fig. 9). Illustrations of the 
mechanisms of the formation of the 1st spike wings and of the corresponding theoretical lines in 
Fig. 15(a). Boxes (a), (b) and (c) illustrate the lines of Eqs. (104), (125) and (148) respectively: 
the corresponding perturbation parameters are (CO/ = 0.39,/! = 0.0077), {(Of = 0.41,/! = 0.00598) 
and ((Of = 0.43, /; = 0.01009) respectively. Resonance separatrices are drawn in red and cyan. The 

dashed lines show the functions /qss(v/) and 1qI^[w)- The black line in (c) is the trajectory of the 
resonant Hamiltonian system (87) which is tangent to both dashed lines. 
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to CO,,, may still give rise to a large variation of action along the trajectory of the 
Hamiltonian system (87). For the odd/even spikes, the boundaries of the chaotic 
layers in the asymptotic limit ^> ^ are formed in this case by the trajectory of (87) 
which is tangent to the lower/upper GSS curves (for the lower/upper layer) or by the 
lower/upper part of the separatrix of (87) generated by the saddle (for the 

upper/lower layer). The overlap of the layers occurs when these trajectories coincide 
with each other, which may be formulated as the equaUty of H in the corresponding 
tangency and saddle: 



for = 1,3,5,..., 
/?(/,, Ws) = H{ll"\wi"^) for j = 2,4,6, . . . , 
-5,), 7,^ 7(4') + 5,). 



(147) 



Note however that, for moderately small 4>, the tangencies may be relevant both to 
the lower layer and to the upper one (see the Appendix). Indeed, such a case occurs 
for our example with <P = 0.2: see Fig. 14(c). Therefore, the overlap of the layers 
corresponds to the equality of H in the tangencies: 



To the lowest order, Eq. (147) and Eq. (148) read as: 



h = h{(Of) 



COf 



2 In 



'4e^ 



(148) 



(149) 
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Eq.(149) 
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Fig. 15 (Color version may be found in the online version of 1431 as Fig. 10). The 1st (a) and 
2nd (b) spike in hg^.{af): comparison between the results of the numerical simulations (the lower 
boundary of the shaded area) and the theoretical estimates. The estimates are indicated by the 
corresponding equation numbers and are drawn by different types of lines, in particular the dashed 
lines represent the explicit asymptote for the solid line of the same color. 
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Both the line (148) and the asymptotic line (149) well agree with the part of the right 
wing of the 1st spike situated to the right from the fold at (Of ~ 0.42 (Fig. 15(a)). 
The fold corresponds to the change of the mechanisms of the chaotic layers overlap. 

If 4> is moderately small while n > 1, the description of the far wings by the 
numerical lines (104) and (148) may be still quite good but the asymptotic lines 
(106) and (149) cannot pretend to describe the wings quantitatively any more (Fig. 
15(b)). As for the minimum of the spike and the wings in its close vicinity, one 
more mechanism may become relevant for their formation (Figs. 15(b) and 16). It 
may be explained as follows. If « > 1, then qn{E) becomes zero in the close vicinity 
(~ ^>2) of the relevant barrier (the upper or lower barrier, in the case of even or odd 
spikes respectively; cf. Fig. 16). It follows from the equations of motion (98) that 
no trajectory can cross the line / = Iq„={). In the asymptotic limit ^> ^ 0, provided 
h is from the range relevant for the spike minimum, almost the whole GSS curve 
is further from the barrier than the line / = Iq„=Q, and the latter becomes irrelevant. 
But, for a moderately small <P, the line may separate the whole GSS curve from 
the rest of the phase space. Then the resonance separatrix cannot connect to the 
GSS curve even if there is a state on the latter curve with the same value of H as 
on the resonance separatrix. For a given (Of, the connection then requires a higher 
value of h: for such a value, the GSS curve itself crosses the line / = /i,„=o. In the 
relevant range of h, the resonance separatrix passes very close to this line, so that 
the connection is well approximated by the condition that the GSS curve touches 
this line (see the inset in Fig. 16); 

5„ ^fif -£'i,2;_i=o for 7=2,4,6,..., 
5,=£,2,_.=o-4'^ for = 3,5,7,.... (150) 

This mechanism is relevant to the formation of the minimum of the 2nd spike at 
= 0.2, and in the close vicinity of the spike on the left (Fig. 15(b)). 



Fig. 16 (Color version may 
be found in the online version 
of (43l as Fig. 11). Amplitude 
of the 3rd Fourier harmonic 
as a function of action (solid 
red line). The dashed black 
line shows the zero level. Its 
intersection with the solid red 
line is marked by the circle. 
The green line indicates the 
value of action where qj, = 0. 
The inset illustrates the line 
(150) in Fig. 15(b): the GSS 
curve touches the horizontal 
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Finally, let us find explicitly the universal asymptotic shape of the spike in the 
vicinity of its minimum. First, we note that the lowest-order expression (146) for 
the spike between the minimum and the fold can be written as the half-sum of the 
expressions (106) and (149) (which represent the lowest-order approximations for 
the spike to the left of the minimum, and to the right of the fold respectively). Thus, 
all three lines (106), (146) and (149) intersect at a single point. This means that, in 
the asymptotic limit <P ^Q, the fold merges with the minimum: (Of and h in the fold 
asymptotically approach Oj and hs respectively. Thus, though the fold is a generic 
feature of the spikes, it is not of major significance: the spike is formed basically 
from two straight lines. The ratio between their slopes is universal. So, introducing 
a proper scaling, we reduce the spike shape to the universal function (Fig. 17): 



h{A(bf) = = 1 - Vl -4e'-2ziw/ w 

w 1 - 0.5934 % for 4 % < 0, 
h{A 6}f) = /i'™') {Aa)f) = \ + Aa)f for 4 % > 0, 

(151) 

1 - Vl — 4e''"2 
= 1 + 4% w 1+0.2034%, 

r h (Of- ©51 

h=—, Aa)f = -^ 0^0, 

hsQ CO,i - (0,0 

where (Oso and h^o are the lowest-order expressions (136) respectively for the fre- 
quency and amplitude in the spike minimum, (Osi is the expression (139) for the 
frequency in the spike minimum, including the first-order correction, and c is a con- 
stant (135). 

In addition to the left and right wings of the universal shape (the solid lines in 
Fig. 17), we include in (151) the function h''f"'^\A6)f) (the dashed line in Fig. 17): 
its purpose is to show, on one hand, that the fold merges asymptotically with the 
minimum but, on the other hand, that the fold is generic and the slope of the spike 



Fig. 17 (Color version may 
be found in the online version 
of (43j as Fig. 12). The 
universal shape of the spike 
minimum (151) (solid lines). 
The dashed line indicates the 
universal slope of the spike 
in between the minimum and 
the fold, which have merged 
in the universal (asymptotic) 
function (151). 
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between the minimum and the fold has a universal ratio to any of the slopes of the 
major wings. 

Even for a moderately small (p, as in our example, the ratios between the three 
slopes related to the 1st spike in the simulations are reasonably well reproduced by 
those in Eq. (151): cf. Figs. 15(a) and 17. It follows from (151) that the asymptotic 
scaled shape is universal i.e. independent of <P (but still assuming the asymptotic 
limit <P 0), « or any other parameter 

The description of the wings of the spikes near their minima, in particular the 
derivation of the spike universal shape, constitutes the third main result of this sec- 
tion. 



4.3 Generalizations and Applications 

The facilitation of the onset of global chaos between adjacent separatrices has a 
number of possible generalizations and applications. We discuss an application in 
Sec. 5, but first list some of generalizations below. 

1. The spikes in hgdcOf) may occur for an arbitrary Hamiltonian system with two 
or more separatrices. The asymptotic theory can be generalized accordingly. 

2. The absence of pronounced spikes at even harmonics IjoJm is explained by the 
symmetry of the potential (84): the even Fourier harmonics of the coordinate, 
q2j, are equal to zero. For time-periodic perturbation of the dipole type, as in 
Eq. (85), there are no resonances of even order on account of this symmetry 
ifTOl |23l l55l ISTl |52l |4ll . If either the potential is non-symmetric, or the additive 
perturbation of the Hamiltonian is not an odd function of the coordinate, then 
even-order resonances do exist, resulting in the presence of the spikes in hgc{(Of) 
at (OfKi2j(o,„. 

3. There may also be an additional facilitation of the onset of global chaos that 
could reasonably be described as a "secondary" facilitation. Let the frequency 
(Of be close to the frequency ©j of the spike minimum, while the amplitude h 
be ~ hs but still lower than hgc{(Of). Then there are two resonance separatrices 
in the / — i/jf plane that are not connected by chaotic transport (cf. Fig. 11(b)). 
This system possesses the zero-dispersion property. The trajectories of the reso- 
nant Hamiltonian (87) which start in between the separatrices oscillate in / (as 
well as in di///df). The frequency dj of such oscillations along a given trajectory 
depends on the corresponding value of ft analogously to the way in which (O de- 
pends on E for the original Hamiltonian //q: 6>{H) is equal to zero for the values 
of H corresponding to the separatrices (being equal in turn to Hsi and Hsu'- see 
Eq. (103)) while possessing a nearly rectangular shape in between, provided the 
quantity 1^^/ — fisu\ is much smaller than the variation of H within each of the 
resonances, 

\Hsl — Hsu I ^ Hi,ar ~ \Hsl — Hel \ ~ — | , (152) 



A New Approach To The Treatment Of Separatrix Chaos And Its Applications 



59 



where Hgi and //g„ are the values of H at the elHptic point of the lower and 
upper resonance respectively. The maximum of 6}{H) in between Hgi and H„, is 
described by the asymptotic formula: 

(»m« 7-^ t!^ ^-;t- (153) 

^ In (//,„,/!//,,-//,„ I) 

If we additionally perturb the system in such a way that an additional time- 
periodic term of frequency 0/ w dim arises in the resonance Hamiltonian, then 
the chaotic layers associated with the resonance separatrices may be connected 
by chaotic transport even for a rather small amplitude of the additional perturba- 
tion, due to a scenario similar to the one described in this paper 
There may be various types of such additional perturbation fyT\. For example, 
one may add to H (85) one more dipole time-periodic perturbation of mixed 
frequency (i.e. « cOm + ^m)- Alternatively, one may directly perturb the angle of 
the original perturbation by a low-frequency perturbation, i.e. the time-periodic 
term in H (85) is replaced by the term 

-/z^cos(a)/f +Acos(fij/-f)), (OfKiO},„, (3/ w «,„. (154) 

Recently discussed physical problems where a similar situation is relevant are: 
chaotic mixing and transport in a meandering jet flow ll30l and reflection of light 
rays in a corrugated waveguide 122] . 
4. If the time-periodic perturbation is multiplicative rather than additive, the reso- 
nances become parametric (cf. [|2ll). Parametric resonance is more complicated 
and much less studied than nonlinear resonance. Nevertheless, the main mecha- 
nism for the onset of global chaos remains the same, namely the combination of 
the reconnection between resonances of the same order and of their overlap in 
energy with the chaotic layers associated with the barriers. At the same time, the 
frequencies of the main spikes in hgc{(Of) may change (though still being related 
to (0,n). We consider below an example when the periodically driven parameter 
i@ in (84). The Hamiltonian is 



i/ = pV2 + (^-sin(^))V2, 

<P = <P[) + hco5{(Oft), 4>o = const<l. (155) 

The term (0 — sin(^))^/2 in H (155) may be rewritten as (^>o — sin(^))^/2 + 
(00 ~ sin(^))/jcos(a)/'f) +/!^cos^(a)/f)/2. The last term in the latter expression 
does not affect the equations of motion. Thus, we end up with an additive per- 
turbation (0() — ?,m{q))hcos{(£)ft). In the asymptotic limit ^>o 0, the nth-order 
Fourier component of the function (^o ^ sin(q') ) can be shown to differ from zero 
only for the orders n = 2,6,10,... Therefore one may expect the main spikes in 



In the case of a 2D electron gas in a magnetic superlattice, this may correspond e.g. to the 
time-periodic electric force applied perpendicular to the direction of the periodic magnetic field 
(49] 1^. 
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Fig. 18 Diagram analogous 
to that in Fig. 9, but for the 
system (155) (with <Po = 0.2). 
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hgc{cOf) to be at frequencies twice larger than those for the dipole perturbation 
(85): 

co^j) « 2a)P « 2(2y- 1)0),,,, = 1,2,3, ... (156) 
This agrees well with the results of simulations (Fig. 18). 

Moreover, the asymptotic theory for the dipole perturbation may immediately 
be generalized to the present case: it is necessary only to replace the Fourier 
component of the coordinate q by the Fourier component of the function {<Po~ 
sin(^)): 

{00 - Mq))n - '1, "TJmX j ^ ^0 ^ (157) 

(cf. Eq. (96) for q„). We obtain: 



21n(^) 8i„(^) 
« = 2,6,10,..., ^0^0, 
where c is given in Eqs. (134) and (135). 

For 00 = 0.2, Eq. (158) gives, for the 1st spike, values differing from the sim- 
ulation data by about 3% in frequency and by about 10% in amplitude. Thus, 
the lowest-order formulae accurately describe the 1st spike even for a moderately 
small <P. 



A New Approach To The Treatment Of Separatrix Chaos And Its Applications 



61 



5. One more generalization relates to multi-dimensional Hamiltonmn systems with 
two or more saddles with different energies: the perturbation may not necessarily 
be time-periodic, in this case. The detailed analysis has not yet been done. 

The paper ||43l presents a rather detailed discussion of possible applications to 
the telectron gas in a magnetic superlattice, a spinning pendulum, cold atoms in an 
optical lattice as well as to problems of noise-induced escape and the stochastic web 
formation. We review briefly in the next section the further development of the latter 
application. 



5 Enlargement of a low-dimensional stochastic web 

The stochastic web concept dates back to the 1960s when Arnold showed lH that, 
in non-degenerate Hamiltonian systems of dimension exceeding 2, resonance lines 
necessarily intersect, forming an infinite-sized web in the Poincare section. It pro- 
vides in turn for a slow chaotic (sometimes called "stochastic") diffusion for infinite 
distances in relevant dynamical variables. 

It was discovered towards the end of 1980s 1541 l7l [8l l9l that, in degenerate or 
nearly-degenerate systems, a stochastic web may arise even if the dimension is 3/2. 
One of the archetypal examples of such a low-dimensional stochastic web arises 
in the ID harmonic oscillator perturbed by a weak traveling wave the frequency of 
which coincides with a multiple of the natural frequency of the oscillator lISTl lSl lSSl . 
Perturbation plays a dual role: on the one hand, it gives rise to a slow dynamics char- 
acterized by an auxiliary Hamiltonian that possesses an infinite web-like separatrix; 
on the other hand, the perturbation destroys this self-generated separatrix, replacing 
it by a thin chaotic layer. Such a low-dimensional stochastic web may be relevant to 
a variety of physical systems and plays an important role in corresponding transport 
phenomena: see ifSTl [8] |55] for reviews on relevant classical systems. In addition, 
there are quantum systems in which the dynamics of transport reduces to that in the 
classical model described above. The latter concerns e.g. nanometre-scale semicon- 
ductor superlattices with an applied voltage and magnetic field lfT5][T6l . 

One might assume that, like the Arnold web, the low-dimensional stochastic web 
described above should be infinite, so that it can provide for transport between the 
centre of the web and states situated arbitrarily far away in coordinate and momen- 
tum. However the numerical integration of the equations of motion shows that this 
is not so: even for a rather non-weak perturbation, the real web is limited to the re- 
gion within a few inner loops of the infinite web-like resonant separatrix (Fig. 19(a)) 
while chaotic layers associated with outer loops are distinctly separated from each 
other and from the web-like chaotic layer formed by the few inner loops. The reason 
is apparently as follows. The single infinite web-like separatrix is possessed by the 
resonant Hamiltonian only in the first-order approximation of the averaging method 
||5l whereas, with the account taken of the next-order approximations, the separa- 
trix apparently splits into many separate complex loops successively embedded into 
each other. Non-resonant terms of the perturbation dress the separatrices by expo- 
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nentially narrow chaotic layers. If the perturbation is not small, the chaotic layers 
manage to connect neighbouring separatrix loops situated close to the centre. How- 
ever, the width of the chaotic layer decreases exponentially sharply as the distance 
from the centre grows i51, 8, 55 1. As a result, the merger between chaotic layers 
associated with neighbouring loops takes place only within the few loops closest to 
the centre, provided that the perturbation is not exponentially strong. 

If the resonance between the perturbation and the oscillator is inexact, or if the os- 
cillator is nonlinear, the splitting between the neighbouring loops is typically much 
larger: it appears even in the first-order approximation of the averaging method 
lISTl |9] \55l . So the number of loops connected to the centre by chaotic transport 
is even smaller 1511 191 1551 than in the case of the exact resonance. 

A natural question arises: how can the perturbation be modified in order for the 
transport to be unlimited or, at least, significantly extended? One of the answers was 
obtained in the very beginning of studies of the low-dimensional webs 1541 171: if the 
perturbation consists of repeated in time short kicks that are also periodic space, and 
if the frequency of the kicks is equal to a multiple of the natural frequency, then a 
so-called uniform web covering the whole of phase space is formed. However such 
a perturbation is absent in many cases and, even where present, the chaotic transport 
is still exponentially slow if the perturbation is weak lISTllSSl . 

It is reasonable then to pose the following question: is it possible to obtain a web 
of form similar to the original one |^ but substantially extended in phase space? 
A positive answer was suggested in l43l and explicitly realized recently 1461 using 
the following simple idea. The chaotic layer in the webs is exponentially narrow 
since the frequency of the non-resonant perturbation of the resonant Hamiltonian is 
necessarily much higher than the frequency of small eigenoscillation in the cell of 
the web-like separatrix ISTl l54l 171 ISl l9l l55l . So we need to modify the perturbation 
in such a way that the resonant Hamiltonian does not change while its perturbation 
contains, in addition to the conventional terms, a low-frequency one. One may do 
this modulating the wave angle with a low frequency or adding one more wave 
with the frequency slightly shifted from the original one. The latter option, together 
with a generalization for the uniform web leading to a huge enhancement of the 
chaotic transport through it, have not yet been considered in detail while the work 
l46ll and the present section concentrate on the former option since it may have 
immediate applications to nanometre-scale semiconductor superlattices in electric 
and magnetic fields ITSl IT6l . 



5.1 Slow modulation of the wave angle 



Fig. 19 demonstrates the validity of our idea. We integrate the equation 

+ ^ = 0. 1 sin[15^ - 4f - /isin(0.02f )] , (159) 
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first for li ~ (i.e. for the conventional case with parameters as in IfSTl [8] |55l ). 
and secondly for /i = 0.1. Although the modulation in the latter case is weak (its 
amplitude is about 63 times smaller than the 2n period of the wave angle which is a 
characteristic scale in this problem), the resultant increase in the size of the web in 
coordinate and momentum is large: by a factor of ~6. 

An analytic theory can be developed to account for these results. It can be gener- 
alized for the off -resonant case 1511 191 155]| too, using the general method developed 
in P3ll44l[38l[39l and described above in the previous sections. 

It is anticipated that the method can also be generalized for uniform webs lISTl 
l54l l55l too, leading to an exponentially strong enhancement of chaotic transport 
through them. 



5.2 Application to semiconductor superlattices 

The works ifTSl [T6ll consider quantum electron transport in ID semiconductor su- 
perlattices (SLs) on the nanometre scale, subject to a constant electric field along 
the SL axis and to a constant magnetic field. The spatial periodicity with a period of 
the nanometre scale gives rise to the onset of minibands for electrons. In the tight- 
binding approximation, the electron motion in the lowest mini-band is described by 
the following dispersion relation for the electron energy E versus momentum p; 
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Fig. 19 The Poincare section for a trajectory of the system (159) with initial state 17 = 0.1, 4 = 
(at instants t„ = nT where T = 2;r/0.02 is the period of the modulation and n — 1,2,3, ...600000) 
for /i = (left panel) and /; = 0. 1 (right panel). A sympletic integration scheme of the fourth order 
is used, with an integration step f,„, = -^^q ~ 1.57 x 10^*, so that the inaccuracy at each step is 
of the order of f?j, Ri x 10^". The left panel coiTcsponds to the example of the conventional case 
considered in I51| [8l [55l . The right panel demonstrates that the modulation, although weak, greatly 
enlarges the web sizes. 
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where x is the direction along the SL axis, A is the miniband width, d is the SL 
period, m* is the electron effective mass for the motion in the transverse (i.e. y — z) 
direction. 

Thus, the quasi-classical motion of electron in an electric field F and a magnetic 
field B is described by the following equation: 

^ = -e{F+[Vp£(p)xB]}. (161) 

where e is the electron charge 

It was shown in flSl that, with a constant electric field along the SL axis 
F — (— Fo,0,0) and a constant magnetic field with a given angle d to the axis 
B = (ficos(0),O,fisin(0)), the dynamics of the z-component of momentum pz re- 
duces to the equation of motion of an auxiliary harmonic oscillator in a plane wave. 
At certain values of the parameters, the ratio of the wave and oscillator frequencies 
takes integer values (like in Eq. (159) with h = 0) which gives rise to the onset of the 
stochastic web, leading in turn to a delocalization of the electron in the x-direction 
and, as a result, to an increase of the dc-conductivity along the SL axis. The experi- 
ment lfT6l appears to provide evidence in favor of this exciting hypothethis. 

At the same time, the finite size of the web and, yet more so, the exponentially 
fast decrease in the transport rate as the distance from the centre of the web in- 
creases, seems to put strong hmitations on the use of the effect. We suggest a simple 
and efficient way to overcome these limitations. Indeed, one can show that, if we 
add to the original (constant) electric field Fq a small time-periodic (ac) component 
FacSin(i2acf)> then the wave angle in the equation of motion of is modulated by 
the term (cf. Eq. (159)): 

/!sm(X2f) = — — sm — H, Q.q=——. (162) 

This allows us to increase drastically the size of the web and the rate of chaotic 
transport through it. For example, for the case shown in Fig. 19, where we have an 
increase of the web size by a factor of 6 x , it is sufficient to add an ac component of 
the electric field with the frequency 0.02 • i2o and an amplitude F^c = 0.1 • 0.02 • Fq 
i.e. an amplitude smaller than that of the original constant field Fq by a factor of 
500 X. 



5.3 Discussion 



We have presented above just initial results on the subject P31 . There are still many 
unsolved interesting problems - 
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1. It can be shown that, in the off -resonance case, there may be a facilitation of 
the onset of global chaos similar to that described in Sec. 4 above, i.e. the crit- 
ical value of the modulation amplitude li required for the onset of global chaos 
between neighbouring separatrix loops possesses deep spikes (minima) as a func- 
tion of the modulation frequency Qac- The detailed theory of this faciUtation has 
yet to be developed. 

2. Our conjecture that, in the resonant case, taking account of the next-order ap- 
proximations of the averaging method could explain the split between different 
separatrix loops, should be proved rigorously. If the corresponding theory is de- 
veloped, it will provide the possibility of calculating both the optimal modulation 
frequency, i.e. that at which the web sizes are maximal, for a given amplitude of 
modulation, and the maximum sizes themselves. 

3. It would be interesting to study the case with an additive perturbation (rather than 
an angular modulation) in detail, both numerically and theoretically. 



6 Conclusions 

We have reviewed the recently developed method for the theoretical treatment of 
separatrix chaos in regimes when it involves resonance dynamics. It has been ap- 
plied both to single-separatrix chaotic layers and to the onset of global chaos be- 
tween two close separatrices. The method is based on a matching between the dis- 
crete chaotic dynamics of the separatrix map and the continuous regular dynamics 
of the resonance Hamiltonian. For single-separatrix chaos, the method has allowed: 

1. Development of the first asymptotic (i.e. for h 0) description of the high peaks 
in the width of the separatrix chaotic layer as a function of the perturbation fre- 
quency, thus describing its dominant feature and, in particular, its maxima. 

2. Classification of all systems into two types, based on the asymptotic dependence 
of the maximum width on the perturbation amplitude h: the maximum width is 
proportional to /iln(l /h) or li for systems of type I or type II respectively. 

For systems with two or more separatrices, the method has allowed us to develop an 
accurate asymptotic theory of the facilitation of the onset of global chaos between 
neighbouring separatrices which occurs at frequencies close to multiples of a local 
maximum in the eigenoscillation frequency as a function of the energy: the local 
maximum necessarily exists in the range between the separatrices. 

Finally, for an oscillator perturbed by a plane wave of frequency equal to or 
close to the frequency of a small eigenoscillation, the method has allowed us to 
suggest how to enlarge substantially the size of the stochastic web using a rather 
weak perturbation, and it promises to provide an accurate theoretical description of 
the enlargement. 
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7 Appendix 

This appendix follows the appendix of the paper 143)1 . The chaotic layers of the 
system (85) associated with the separatrices of the unperturbed system (84) are de- 
scribed here by means of the separatrix map. To derive the map, we follow the 
method described in lf55l . but the analysis of the map significantly differs from for- 
merly existing ones ll23l l55l ISTl l52l l29l (cf. also the recently published paper 1351 
where the analysis of the map has some similarity to ours but still differs signifi- 
cantly). Using our approach, we are able to calculate the chaotic layer boundaries 
in the phase space (rather than only in energy), throughout the resonance frequency 
ranges, and we can quantitatively describe the transport within the layer in a manner 
different from existing ones (cf. l29l[32ll and references therein). 



7.1 Lower chaotic layer 



We now present a detailed consideration of the lower chaotic layer. The upper 
layer may be considered in a very similar way. 

7.1.1 Separatrix map 

A typical form of trajectory ^(f) close to the inner separatrix (that corresponding 
to the lower potential barrier) is shown in Fig. 20. One can resolve pulses in ^(f ). 
Each of them consists of two approximately antisymmetric spike^. The pulses are 
separated by intervals during which |^| is relatively small. In general, successive 
intervals differ between each other Let us introduce the pair of variables E and (p: 

E = Hq, (p = (Oft + (Pa, (163) 

where the constant % may be chosen arbitrarily. 



Spikes correspond to motion over any of the minima of the potential, first in one and then (after 
the reflection from one of the upper barriers) in the opposite direction. If <P is small, then the spikes 
within the pulse are separated by long intervals since the reflection point is situated close to the top 
of the upper barrier, where the motion is slow. 
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Fig. 20 Schematic example 
of the time dependence of 
the velocity of the perturbed 
system (thick solid line) in 
the case when the energy of 
motion varies in the close 
vicinity of the top of the lower 
potential barrier The dashed 
line marks the zero level of the 
velocity. Pulses of the veloc- 
ity are schematically singled 
out by the parallelograms 
(drawn by a thin solid line). 
The two sequences of time 
instants (...,fi_i,f;,f,+i, ...) 
and (...,?;_, cor- 
respond to beginnings and 
centers of the pulses, respec- 
tively. 




ti-i t'i-i 



The energy E changes only during the pulses of q[t) and remains nearly un- 
changed during the intervals between the pulses, when |^(f)| is small 1551 . We as- 
sign numbers / to the pulses and introduce the sequences of (£, , (j?,) corresponding 
to the initial instants f,- of the pulses. In such a way, we obtain the following map (cf. 
1551): 



£,■+1 -E, + AEi, <p,+i - <p,- + ^^^^ , 

AEi = h I dt q{t)cos{cOft), (164) 
Ji th pulse 

where jjj puj^g means integration over the ith pulse. Before deriving a more explicit 
expression for AEi, we make two remarks. 

1. Let us denote with the instant within the ith pulse when q is equal to zero 
(Fig. 20). The function q{t — f ■) is an odd function within the /th pulse and it is 
convenient to transform the cosine in the integrand in AEi (164) as 



cos(aj/f) = cos(a)/(f-f,0 + fiJ/fO = cos(tO/-(f-f,'))cos(a)/f')-sin(a)/(f-f-))sin(a)/f/), 

and to put % = (Of{t\ — t^), so that (p; = cOft'-. 
2. Each pulse of q contains one positive and one negative spike. The first spike can 
be either positive or negative. If E changes during the given nth pulse so that its 
value at the end of the pulse is smaller than then the first spikes of the ith 
and (i + l)st pulses have the same signs. On the contrary, if E at the end of the 
;th pulse is larger than fi^'', then the first spikes of the /th and {i+ l)st pulses 
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have opposite signs. Note that Fig. 20 corresponds to the case when the energy 
remains above /s^'' during the whole interval shown in the figure. This obviously 
affects the sign of AEi, and it may be explicitly accounted for in the map if we 
introduce a new discrete variable d; = ±1 which characterizes the sign of q at 
the beginning of a given ith pulse, 

a; = sign(^(r,)) , (165) 
and changes from pulse to pulse as 

CT,+i =c7,sign(4''-£,+ i) . (166) 
With account taken of the above remarks, we can rewrite the map (164) as follows: 



£,+ 1 = Ei + c7,/ie('""') sin((j!),), (167) 
a)/7r(3-sign(£,-+i-4'^)) 

(Pi+l = (Pi + — TT—T^ ^ , 

Oi+i = (J,- sign (£'^'' -£■,■+!), 

£(/ow) ^ = -cj,- / At q{t - 1[) sin(a)/(f - 1[)) 

J ith pulse 

f'i+l 

K-20i / dr^(f-f-)sin((B/(f-r,')). 
Jt' 



A map similar to (167) was introduced in ||53l , and it is often called the Zaslavsky 
separatrix map. It was re-derived mathematically rigorously in ISTl : see also the 
recent mathematical review ||29l . The latter review also describes generalizations of 
the Zaslavsky map as well as other types of separatrix map. The analysis presented 
below relates immediately to the Zaslavsky map but it is hoped that it will prove 
possible to generalize it for other types of separatrix maps too. 

The variable gC"") introduced in (167) will be convenient for the further calcu- 
lations since it does not depend on ; in the lowest-order approximation. A quantity 
like 5/ = /i I e is sometimes called the separatrix split fSTl since it is convention- 
ally assumed that the maximal deviation of energy on the chaotic trajectory from the 
separatrix energy is of the order of 5/ 123] |55] |5l] |52l . As in the main text, we shall 
use this term, but we emphasize that the maximal deviation may be much larger 

In the adiabatic limit (Of 0, the excess of the upper boundary Zs^''' of the lower 

layer over the lower barrier e'^P does not depend on angle and is equal to 2nh (cf. 
lfT4l ). But (Of relevant for the spike of hgc{(Of) cannot be considered as an adiabatic 
frequency, despite its smallness, because it is close to (£),„ or to its multiple while all 
energies at the boundary lie in the range where the eigenfrequency is also close to 

cOf^{2j-\)co,„^{2j-\)co{E^^), = 1,2,3,... (168) 
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The validity of (168) (confirmed by the results) is crucial for the description of the 
layer boundary in the relevant case. 

7.1.2 Separatrix split 

Let us evaluate g ('""') explicitly. Given that the energy is close to fi^' , the velocity 
q{t — ?•) in (167) may be replaced by the corresponding velocity along the 

separatrix associated with the lower barrier, qs"^''\t — ?■), while the upper limit of 
the integral may be replaced by infinity. In the asymptotic limit <P -^0, the interval 
between spikes within the pulse becomes infinitely long^^ so that only the short 
(~ Wq"') intervals corresponding to the spikes contribute to the integral in gC'™') 
(167). In the scale COj^, they may be considered just as the two instants: 

In the definition of gC'™) (167), we substitute the argument of the sine by the corre- 
sponding constants for the positive and negative spikes respectively: 



e('-)«2sinf^') / df^f"')(f-r;)«2;rsinfj^y (170) 

\4-C0„, J Jpositive spike \4C0,„ J 

0^0. 

In the derivation of the first equality in (170), we have also taken into account that 
the function qs''^''\x) is odd. In the derivation of the second equality in (170), we 
have taken into account that the right turning point of the relevant separatrix is the 
top of the lower barrier and the distance between this point and the left turning point 
of the separatrix approaches n in the limit 0—^0. 

For the frequencies relevant to the minima of the spikes of hgc{cOf), i.e. for O/ = 

Co]-'' w (2 j — 1 )ft),n, we obtain: 



^ 2;rsin ((27- 1)|) = V27r(-1)[^], 
1,2,3,..., cp^O. (171) 

For moderately small 0, it is better to use the more accurate formula: 

£('™')(a)^) = 2 Tdf 4S'""'(f)sin(a)/f), (172) 



where the instant t — corresponds to the turning point of the separatrix to the left of 
the lower barrier, i.e. qi'"*^ (f = 0) = while qi'"" ^ > for all f > 0. The dependence 
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by Eq. (172) is shown for <P ~ 0.2 in Fig. 21(a). For small frequencies, 



the asymptotic formula (170) fits well the formula (172). 



7.1.3 Dynamics of the map 

Consider the dynamics of the map (167) when (Of is close to the spikes' minima: 
(Of « ncOrn where « = 2y — 1 while j = 1,2,3, .. .. Let the energy at the step / = — 1 
be equal to Ej^^\ The trajectory passing through the state with this energy is chaotic 

since {a){E))^^ diverges as £ — > ii^'^ and, therefore, the angle (p^i is not correlated 
with the angle on the previous step (p_2- The quantity (7_i is not correlated with 
a-2 either. Thus, sm{(p_i) may take any value in the range [—1,1] and (7_i may 
equally take the values 1 or -1. Therefore, the energy may change on the next step 
by an arbitrary value in the interval [~'h\e^'"^'^\,h\e^'""''' \]. Thus, £0 — may have 
a positive valuj^ '--^ /z|e''™^ |. Then, the approximate equality nco{Eo) ~ ft),,, holds, 
provided that the value of li is from the relevant range. Allowing for this and recall- 
ing that we are interested only in those realizations of the map such that £0 > 
the relevant realization of the map i = — I — > / = may be written as: 

£0 = 4'^ + c^-i/ie*'""' sin((p_i) = sin((p_i)|, 
(Po ~ (P-i +nn, 

(Jo =-c7-i. (173) 



7 




7 




Fig. 21 Theoretical estimates for the nonnahzed separatrix spUt (for <P = 0.2) as a function of 
the perturbation frequency, for the lower and upper layers in (a) and (b) respectively. The solid 
lines are calculated from Eqs. (172) and (204) for (a) and (b) respectively, while the dashed lines 
represent the asymptotic expressions (170) and (205) respectively. 



The latter is valid for any (p_ 1 except in the close vicinity of multiples of Tt while the state £0 , (jOo 
(167) in the latter range of turns out irrelevant to the boundary, as shown further down. 
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One may expect that further evolution of the map will, for some time, approxi- 
mately follow the trajectory of the system (87) with the initial energy £0 (173) and 
an arbitrary (p_i and initial slow angle iff somehow related to (po ~ <P-i + n^- Let us 
prove this explicitly. 

Consider two subsequent iterations of the map (167): 2/ ^ 2/ + 1 and 2/ + 1 ^ 
2i + 2 with an arbitrary / > 0. While doing this, we shall assume the validity of 
(168) (it will be clarified below when this is true) from which it follows that: (i) 
(o{Eii^i) w co{Ei^), (ii) (pk+i ~ (pi^ K iiK = (2j — l);r. It will follow from the results 
that the neglected corrections are small in comparison with the characteristic scales 
of the variation of E and <p (cf. the conventional treatment of the nonlinear resonance 
dynamics QUI |23l |55l Ell |52l gl]). Furthermore it follows from (167) that, while 
the energy remains above the barrier energy, Oj^ oscillates, so that 021 = do and 
CT2i+i = -c7o- Then, 



Eii+i =E2i + aohe'^'"'"hm{(p2i) 

(Of 



(P2i+l = (P2i 



(0{E2i+l) 



n ~ (p2i + nn + n 



(Of-n(o{E2i) 

(0{E2i) ' 



(174) 



E2i+2 = E2i+i - OQhe^'""^ M92,+\) = 

= £2,-+! + CTohe^''"'^ sin((p2,'+i - nn) w £2,- + CJQlhe^''"'^ sin(^2/), 

COf (Of -n(0(E2i) 

(p2i+2 = (p2i+i + w (p2i + 2Kn + In ' \ (175) 

(0{E2i+2) 0}{E2ij 

(the second equality in the map for £21+2 takes into account that n is odd so that 
sin(^ — nn) = — sin(^).) 

The quantity (p2i+2 — (P2i — 27r« is small, so the map 2/ ^ 2/ + 2 (175) may be 
approximated by differential equations for £2, and 92; = (Pii — 2nni: 



j^.M."».sin(fe), -^^.-|_^,,-,..^E,,)). (176, 

= (P2i-27zni. 

Let us (i) use for gC'™') the asymptotic formula (171), (ii) take into account that 
the increase of / by 1 corresponds to an increase of time by n/(o{E), and (iii) trans- 
form from the variables {E,(p) to the variables (/,v/ = «;r(l — (7o)/2 — ^). Equations 
(176) reduce then to: 

^ = -/7\/2(-l)[3]sin(v/), ^=n(0-(Of, (177) 
1 - Co 

V/ = «7r— (p, n = lj-\. 
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Equations (177) are identical to the equations of motion of the system (87) in 
their lowest-order approximation, i.e. to equations (98) where q„ is replaced by its 
asymptotic value (96) and the last term in the right-hand part of the second equation 
is neglected, being of higher order in comparison with the term n(0 — (Of. 

Apart from the formal identity of Eqs. (177) and (98), in (177) and x^f in (98) 
are identical to. Necessarily t[ corresponds to a turning point (see Fig. 20) while the 
corresponding is equal to Ini or tt + 2ni for the right and left turning points 
respectively (see (87)) i.e. = Ini + ;r(l — c7,)/2, so that ij/jgg) = — O/f = 
«7r(l-(j)/2-^ = VA(i77). 

The relevant initial conditions for (177) follow from (173) and from the relation- 
ship between ^ and (p: 

7(0) = I{E = + /^^^l sin(v/(0))|), (178) 
while Vi^(O) = nn{\ — ctq) /2 — <Po may be an arbitrary angle from the ranges where 

(-l)l"/41sin(v/(0)) <0. (179) 

For moderately small <P, it is better to use the more accurate dynamic equations 
(98) instead of (177) and the more accurate initial value of action instead of (178): 

/(0)=/(£=4'^ + 5,|sin(v/(0))|), di^h\e'^'°% (180) 

with eC™') calculated by (172). 

We name the quantity 5/| sin(ij/)| the generalized separatrix split (GSS) for the 
lower layer. Unlike the conventional separatrix spUt 5/ ifSH , it is angle-dependent. 
The curve /(i//) = I{E = e\^^ + 5/|sin(\J/)|) may be called then the GSS curve for 

the lower barrier and denoted as IqIs{w)- 

Finally, let us investigate an important issue: whether the transformation from 
the discrete separatrix map (i.e. (174) and (A14)) to the differential equations (176) 
is valid for the very first step and, if it is so, for how long it is valid after that. The 
transformation is valid as long as (o{Eii) w ncOf i.e. as long as E/^ is not too close to 

the barrier energy E^^\ At the step k = 0, the system stays at the GSS curve, with 
a given (random) angle \ff{Q) from the range (179). Thus, at this stage, the relation 
(168) is certainly valid (for the relevant range of h and for any angle except for the 
close vicinity of the multiples of n). The change of energy at the next step is positive 
too: 

£1 -£0 = c7o/^e''™'''sin(^o) ~ -C7_i/!e*'""'^ sin(9_i -nn) = 
^ (7_i/i£(''™' sin(^_i) =£0 -£-1 > 0. 

This may also be interpreted as a consequence of the first equation in (177) and of 
the inequality (179). 

Hence, (168) is valid at the step ^ = 1 too. Similarly, one can show that E2—E1 > 
0, etc. Thus, the transformation (174, 175)^(176) is valid at this initial stage indeed. 
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and the evolution of (E, (p) does reduce to the resonant trajectory (14) with an initial 
angle from the range (179) and the initial action (180). This lasts until the resonant 
trajectory meets the GSS curve in the adjacent n range of iff i.e. at t such that the 
state (/(f ), V/(f )) satisfies the conditions: 

n0-4ss(r(0), [v/(f)/^]-[V>(0)/?r] = l. (181) 

At this instant, the absolute value of the change of energy E/t in the separatrix map 
(174) is equal to £4. — E^^^^ (just because the state belongs to the GSS curve) but 
the sign of this change is negative because the sign of sin(^^) is opposite to that of 
sin(^). Therefore, at the step ^+1, the system gets to the separatrix itself, and the 
regular-like evolution stops: at the next step of the map, the system may either again 
get to the GSS curve but with a new (random) angle from the range (179), and start 
a new regular-like evolution as described above; or it may get to the similar GSS 
curve below the barrier and start an analogous regular-like evolution in the energy 
range below the barrier, until it stops in the same manner as described above, etc. 

This approach makes it possible to describe all features of the transport within 
the chaotic layer In the present context, it is most important to describe the upper 
outer boundary of the layer. 

7.1.4 Boundary of the layer 

We may now analyze the evolution of the boundary of the layer as h grows. Some 
of the stages of the evolution are illustrated by Figs. 13, 14 and 22. 

It follows from the analysis carried out in the previous subsection that any state 
(in the / — ij/ plane) lying beyond the GSS curve but belonging to any trajectory fol- 
lowing the equations (98) which possesses common points with the GSS curve be- 
longs to the chaotic layer: the system starting from such a state will, sooner or later, 
reach the separatrix where the chaotization will necessarily occur Therefore, the up- 
per boundary of the chaotic layer coincides with the trajectory following equations 
(98) with the initial action (180) and an initial angle 1/^(0) from the range (179) such 
that the trajectory deviates from the barrier energy by more than does a trajectory 
(98)-(179)-(180) with any other initial angle. There may be only two topologically 
different options for such a trajectory: either it is tangent to the GSS curve, or it is 
the separatrix trajectory which intersects the GSS curve (some schematic examples 
are shown in Figs. 22(a) and 22(b) respectively; some real calculations are shown in 
Figs. 13 and 14). 

1 . Relatively small h 

Consider first values of h which are large enough for the condition (168) to be sat- 
isfied (the explicit criterion will be given in (192)) but which are smaller than the 
value he) = ha- {(Of) determined by Eq. (125) (its meaning is explained below). The 
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Fig. 22 (Color version may be found in tlie online version of t43| as Fig. 16). A schematic figure 
illustrating the formation of the boundary of the lower chaotic layer for h < hi!,- {a)f) in the ranges 
of (Of relevant to (a) odd, and (b) even spikes. The dashed magenta line shows the GSS curve in 

the energy-angle plane: E{iff) = = + 5;| sin(i/ir)|. Green lines show examples of those 

trajectories (98) which have points in common with the GSS curve. One of them (shown by the 
thick green line) relates to the formation of the upper boundary of the lower chaotic layer: in (a), 
the boundary is the trajectory tangent to the GSS curve; in (b), the boundary is the upper part of 
the separatrix generated by the saddle ".v". Yellow dots indicate the relevant common points of the 
GSS curve and the thick green line. They have angles ±1//,''' and energy in the case (a), and 
angles ±1//,-'' and energy s''' in the case (b). The minimum and maximum deviation of energy 

on the boundary from the baiTier energy are denoted as S^^^^ and Smlx respectively. The maximum 
deviation on the GSS curve is equal to 5/. 



further analysis within this range of h differs for the ranges of O/ relevant to odd 
and even spikes, and so we consider them separately. 

A. Odd spikes 

The relevant frequencies are: 

COf^n(0,n, « = 2;-l, 1,3,5,... (182) 

Let us seek the state {//''', '//'/'''} (with within the range ]0, n[) where the res- 
onant trajectory is tangent to the GSS curve. With this aim, we equalise both the ac- 
tions and the derivatives of both curves. The equality of actions immediately yields 
//'' via xfrj'^: = I{E = = Iclsift'^)- The derivative along the GSS curve 
is obtained by differentiation of /gss(v/). The derivative along a resonant trajectory 
can be found by dividing the first dynamic equation in (14) by the second one. Sub- 
stituting the expression of via i//''' into the equality of the derivatives, we obtain 
a closed equation for and its solution immediately gives us the relevant V/(0): 
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'('-''|cos(r/'') (l-;^-/'^cos(v//")) +,„(£)stn(v//") 
= 0, 

=4''+/i|e*'''^''|sin(v//''), f^'^ e [0,7r], 
« = 2j-l, 1,3,5,..., V/(0) = V//". (183) 

A careful analysis of the phase space structure shows that, in the present case (i.e. 
when h < hc){cOf) while j is odd), there is no separatrix of the resonant Hamiltonian 
(4) which would both intersect the GSS curve and possess points above the tangent 
trajector}0. Thus, for this range of h, the outer boundary of the chaotic layer is 
formed by the trajectory following the dynamical equations (98) with the initial 
angle given by (183) and the initial action by (180) (Fig. 22(a)). 

Let us find the lowest-order solution of Eq. (183). We neglect the term 1 — 
(Of/{nCo{E)) (the result will justify this) and use the lowest-order expression for 
the relevant quantities: namely, Eqs. (171) and (96) for gC"*') and q,, respectively, 
and the lowest-order expression for dq„/dE which can be derived from Eq. (95): 



dq„{E) _ 



dE A.fif J7 



4V2(E-El')ln{0 



« = 2y-l, £-4''<^^0. (184) 



Then Eq. (183) reduces to the following equation 

81n(0-i) 

The lowest-order solution of (185) in the range ]0, 7z[ is 



.l-(-l)[5] 



It follows from the definition E,'^ (183) and from (186) that the lowest-order ex- 

f 17(0 17(1) • 
pression for £, is 

£«-4^) = 5,stn(v//'') = ^-==A^. (187) 

The next step is to find the minimum value of the energy on the boundary of the 
layer, fi^';,,. It follows from the analysis of the dynamical equations (98) that the 
corresponding angle ij/mm is equal to if sign(^2;-i) > (i.e. j — 1,5,9, . . .) or to 
n if sign(g'2;_i) < (i.e. j = 3,7, 1 1, . . .): cf. Fig. 8(a). Given that the Hamiltonian 



^■^ For odd numbers j > 3, there are separatrices which lie in the range of E where ft)(£) <C (0,„ i.e. 
much closer to the barrier than the tangent trajectory: these separatrices ai'e therefore irrelevant. 
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(87) is constant along any trajectory (98) while the boundary coincides with one 

such trajectory, the values of the Hamiltonian (87) in the states {/(^j^jn), W ~ Vmin} 

and {//'' , 1//'' } should be equal to each other. In explicit form, this equality may be 
written as 



d£( 1--^) L,(£/")cos(v//'V(-l)[5]^„(£«J) =0. (188) 



jijC' V 'ME) 

mm ^ ^ ^ 

Let us find the lowest-order solution of Eq. (188). Assume that e'^^^ still be- 
longs to the range of E where Co{E) w (0,n (the result will confirm this assump- 
tion). Then the integrand in (188) goes to zero in the asymptotic limit 0^0. 
Hence the integral may be neglected (again, to be justified by the result). The 
remaining terms in Eq. (188) should be treated very carefully. In particular, it 
is insufficient to use the lowest-order value (96) for q,, since it is the difference 
between q„{EP) and qn{E^^^) that matters. Moreover, the approximate equality 

qn{Ei'^)~qn{El2n) ~ dq„{E^'^) /dE{ \E/ '' ~ eIIIj does not apply here either since, 
as follows from Eq. (184), the derivative d^„(£')/d£' may vary strongly in the range 

[eIIJ,,eP] if -£l!L)/(£mL-4") > 1 (agai"' to be justified by the result). 
That is why it is necessary to use the more accurate expression (95) for q„. Allowing 

for the asymptotic expression (186) of i///'' and keeping only the lowest-order terms, 
one can finally reduce Eq. (188) to the relation 




^ 4- (189) 

Substituting here the asymptotic value of e/'^ (187), we obtain the final lowest-order 
expression for the minimum deviation (along the boundary) of the energy from the 
barrier: 



It is necessary and sufficient that the condition Co{E) « co,,, is satisfied at the minimal 
and maximal energies of the boundary to ensure that the second equality in (168) 
holds true, i.e. that Co{E) is close to co,,, for all points of the boundary. 
At the minimal energy, this condition is 

a),„-co(4'' + 5^1,)««»,. (191) 

Eq. (191) determines the lower limit of the relevant range of h. The asymptotic form 
of (191) is: 
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In I 



T, 



<1. (192) 



ln(l/0) 

We emphasize that any h of the order of h^Q (136) satisfies this condition. In the 
asymptotic limit <P ^0, the left-hand part of Eq. (192) goes to zero. 

As for the maximal energy, it may take values up to the energy of the lower saddle 
"si", i.e. Esi (102). Obviously, (168) is valid at this saddle, too. 

B. Even spikes 

The relevant frequencies are: 

«ee2;-1, = 2,4,6,... (193) 

In this case, q„{E) and dq„{E)/dE have different signs for all E within the rel- 
evant range (i.e. where Co{E) w (0,„, qn{E) ~ ^„(£',„)): cf. (96) and (184). Then, in 
the asymptotic limit <P ^Q, Eq. (183) for the tangency does not have any solution 
for ij//'' in the relevant rang There may only be solutions very close to some 
of n integers, and the corresponding energies E^^ are then very close to ii^'' i.e. 
a)(£/'^) ^ 0),,,: therefore they are irrelevant. 

At the same time, unlike for the odd spikes, there exists a saddle with an angle 

Ws~T^ 2^' ^""^^^ 

while the energy (which may be found as the appropriate solution of Eq. (99)) lies 
in the relevant vicinity of the lower barrier (Fig. 22(b)). In the lowest-order approx- 
imation, this saddle energy is: 

+ 5. = -^--^. (195) 

2\/2ln(ln(4e/^')) 

This saddle (denoted in Fig. 22(b) as "i") generates a separatrix. Its upper 
whiskers go to the similar adjacent saddles (shifted in i// by In). In the asymptotic 
limit 0^0, the upper whiskers are much steeper than the GSS curve and hence 
they do not intersect i@ The lower whiskers do intersect the GSS curve and, more- 
over, two intersections lie in the relevant energy range (Fig. 22(b)). Let us show this 
explicitly. We write the expression for the Hamiltonian (87) in the relevant vicinity 
of the barrier energy (i.e. where a),„ — (i){E) <C Om), keeping in the expression both 
the lowest-order terms and the terms of next order (in particular, we use Eq. (95) for 



In case of a moderately small <t>, tangency may exist in the relevant range of energies. The 
boundary of the layer is then formed by the tangent trajectory. 

In case of a moderately small they may intersect the GSS curve. Then, the tangent trajectory 
lying above the separatrix necessarily exists, so the boundary of the layer is formed by this tangent 
trajectory. 
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q„{E) and take into account that < V2 — nqn{E) <^ \f2 for the relevant range of 



(\\ «5ln(2^) / n% \ 2<P /4e 



21n(f) V 21n(f); TT \<P 

0),,, -»(£■ + 5) <c ft),,,. (196) 

The Hamiltonian H should possess equal values at the saddle "i" and at the in- 
tersections of the separatrix and the GSS curve. Let us denote the angle of the inter- 
section in the range ]0,n[ as xj/p, and let us denote the deviation of its energy E^''^ 
fromii^'' as 5-'^ = 5/ sin(i/if-''). 

Assuming that | xfrj'^ — ij/y'' | <c 1 (the result will confirm this) so that 

cos(v>,<")«(-l)["/^l(l-(r*"-V^^")V2)« 

« (-l)l"/41(i _ (5/'V5,)V2) « - {8P/hy-/4), 

the equality of the values of H is: 



In 



,S,,„|»V«<'.,„fl±)U,.v^-IV^-(^. (,97) 
21n(f)l^' \Sj ■ \sp)) 8 !„(*) 2V2I, 

Let us assume that, in the asymptotic limit ^> ^ 0, S/'' ^ 8s (the result will 
confirm this). Then the left-hand part is asymptotically smaller than the first term 
in the right-hand part. So, Eq. (197) implies, in the asymptotic limit, that the right- 
hand side equals zero. Expressing h via 5j from Eq. (195), we finally obtain a closed 
transcendental equation for 5s /Sj-'^: 

In the asymptotic limit ^> ^ 0, the quantity A diverges and, hence, the lowest- 
order asymptotic solution of Eq. (198) reads as 
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Substituting here the expression (195) for ds and the expression (198) for A, we 
obtain: 



1 /«;rln(ln(^)) 
^i^h-J (200) 



(As 

Thus, we have proved the following asymptotic properties of the separatrix gen- 
erated by the saddle "i": (i) the lower whiskers of the separatrix do intersect the GSS 
curve in the relevant range of E (i.e. where the resonant approximation is valid); and 
(ii) the upper whiskers of the separatrix do not intersect the GSS curve (there is no 
solution of Eq. (197) in the range 5/'' > 5s). The former property confirms the self- 
consistency of the asymptotic theory for even spikes; the latter property means that 
the upper outer boundary of the lower chaotic layer is formed by the upper whiskers 
of the separatrix generated by the saddle "s" . 

Finally, we note explicitly that the minimal (along the boundary) deviation of 
energy from the barrier energy occurs exactly at the saddle "s", i.e. 

C = (201) 



2. Relatively large h. 

As h grows, the boundary of the layer rises while the lower part of the reso- 
nance separatrix, on the contrary, falls. They reconnect at the critical value of h, 
h^cr = he) {(Of), determined by Eq. (125), which may be considered as the absorp- 
tion of the resonance by the chaotic layer. If h grows further, then the GSS curve 
and the resonance separatrix intersect. As a result, the trajectory starting from the 
state of angle (183) and action (180), for odd spikes, or from the saddle 'V, for 
even spikes, is encompassed by the resonance separatrix. So, it no longer forms the 
outer boundary of the layer. Rather it forms the inner boundary i.e. the boundary 
of the main island of stability inside the layer, repeated periodically in iff with a 
period 2n (cf. analogous islands in the upper layer in Fig. 13). Unless the lower 
chaotic layer reconnects with the upper one, the outer boundary of the lower layer 
is formed by the upper part of the resonance separatrix. The relevant initial angle 
1/(0) on the GSS curve corresponds to the intersection of the GSS curve with the 
resonance separatrix (cf. the analogous situation for the upper layer in Fig. 13). 
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7,2 Upper chaotic layer 

The upper chaotic layer may be treated analogousl}0 to the lower layer We present 
here only the results. 

Similarly to the lower-layer case, one may consider the ranges of relatively small 
h (namely, smaller than /i'"^ = hcr{(Of) determined by Eq. (126)) and relatively 

large h (i.e. Ii > he')). In the former range, the formation of the boundary occurs in 
a manner which is, in a sense, opposite to that for the lower-layer case. For even 
spikes, the lower outer boundary is formed by tangency while, for odd spikes, it is 
formed by the lower part of the separatrix generated by the saddle "s", analogous to 
the saddle "i"in the lower-layer casj^ 

So, for even spikes, the angle of tangency i///"' is determined by: 



e("'^)|cos(v/r) ( 1 



^^-/i^^cos(v/; ') ) -^„(£)sm(v// ') 



= 0, 

-/.|e("'^)|sin(#) v^/"' £ [0, ;r] , 
« = 2;-l, ;=2,4,6,..., w{Q) = wi"\ 



E=E, 



(202) 



and 1//, determines the tangency energy: 

£/"^=£f' -/2|e(«'')|sin(vA/''^), 
where the quantity e'"''^ is described by the formula 



(203) 



=2 / df ^r''(f)cos(«/f) 



up). 



where q'"''^ (f ) is the time dependence of the velocity along the separatrix associated 



(204) 



with the upper barrier and the instant f = is chosen so that q'"''^ [t = 0) is equal 
to the coordinate of the lower barrier while qf^^ > for f G [0,°°[. The dependence 



{(Of) 



in Eq. (204) is shown for <P = 0.2 in Fig. 21(b). 



The asymptotic form of Eq. (204) is 



For any AC-driven spatially periodic Hamiltonian system, the upper energy boundary of the 
layer associated with the unbounded separatrix diverges in the adiabatic limit (Of —> |42l . How- 
ever, this divergence is not relevant for the present problem for the following reasons. The lower 
chaotic layer relates to the hounded separatrix while, for the upper (unbounded) layer, it is the lower 
boundary of the layer which is relevant for the onset of global chaos in between the separatrices. 
Moreover, even for the upper boundary of the upper layer, the divergence is not yet manifested for 
the driving parameters (/j, cOf) in the vicinity of the spikes minima (cf. 1421 ). 

This tangency may exist for a moderately small <P. The boundary is then formed by the tangent 
trajectory rather than by the separatrix: see an example in Fig. 14(c). 
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^{up) ^ gM (gj^) ^ 27rcos (205) 

For (Of = 0)]-'^ « (2y — l)co,„, Eq. (204) reduces to 

eM(ejP) ^ 2;rcos ((2j- 1)|) = V2n{-l)\-^\ 
= 1,2,3,..., f^O. (206) 

The lowest-order solution of (202) is given in Eq. (121), so that E,"^ is approximated 
by Eq. (122). The maximal energy on the lower boundary of the layer corresponds 
to v/(f) = TTif 7 =2,6,10,... or if ; = 4,8, 12, . . . and is determined by Eq. (123). 
The asymptotic value of the minimal deviation from the upper barrier of the energy 
at the boundary, dj^^^, is given in Eq. (124). 

For odd spikes, the boundary is formed by the lower part of the separatrix gener- 
ated by the saddle "s ". The angle of the saddle is given in Eq. (117) while the devi- 
ation of its energy from the barrier is approximated in lowest-order by Eq. (118). 

As li grows, the boundary of the layer falls while the upper part of the upper 
resonance separatrix rises. They reconnect at li = h["^ = /i'"' (ft)/-), as determined by 
Eq. (126), which may be considered as the absorption of the resonance by the layer 

For larger h, the boundary of the layer is formed by the lower part of the upper 
resonance separatrix (Fig. 13), unless the latter intersects the lower GSS curve (in 
which case, /i'"' marks the onset of global chaos). 
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